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Abstract 

In our companion paper we identified a complete set of manifestly gauge-invariant observables for general 
relativity. This was possible by coupling the system of gravity and matter to pressureless dust which plays the 
role of a dynamically coupled observer. The evolution of those observables is governed by a physical Hamiltonian 
and we derived the corresponding equations of motion. Linear perturbation theory of those equations of motion 
around a general exact solution in terms of manifestly gauge invariant perturbations was then developed. 

In this paper we specialise our previous results to an FRW background which is also a solution of our 
modified equations of motion. We then compare the resulting equations with those derived in standard cos- 
mological perturbation theory (SCPT). We exhibit the precise relation between our manifestly gauge-invariant 
perturbations and the linearly gauge-invariant variables in SCPT. We find that our equations of motion can be 
cast into SCPT form plus corrections. These corrections are the trace that the dust leaves on the system in 
terms of a conserved energy momentum current density. It turns out that these corrections decay, in fact, in 
the late universe they are negligible whatever the value of the conserved current. 

We conclude that the addition of dust which serves as a test observer medium, while implying modifications 
of Einstein's equations without dust, leads to acceptable agreement with known results, while having the 
advantage that one now talks about manifestly gauge-invariant, that is measurable, quantities, which can be 
used even in perturbation theory at higher orders. 
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1 Introduction 



In our companion paper [T] we combined the framework of relational observables [21 [3] with the Brown 
- Kuchaf mechanism [3] in order to cast general relativity (including all known physical matter) into the 
form of an ordinary Hamiltonian system, with a true gauge invariant Hamiltonian which only depends on 
the gauge invariant geometry and matter degrees of freedom. From that foundation we then developed a 
manifestly gauge-invariant perturbation theory. 

The motivation [5] was the trivial observation that general relativity is a gauge theory (the gauge 
group being the diffeomorphism group). Therefore, neither metric nor matter fields are directly observable 
since they are not gauge invariant and likewise Einstein's equations do not describe physical evolution of 
observables but rather the behaviour of non - observables under gauge transformations. On the other hand, 
we would like to think of Einstein's equations as describing the physical time evolution of observables. The 
latter point of view is often justified by the argument that the universe could be filled, in principle, by 
geodesic test observers who serve as a material reference system. The achievement of [1] is to implement 
that idea field theoretically by adding pressureless dust to the system which serves as the medium of test 
observers. What we did in [1] is to extend that idea further and thus to arrive at a fully gauge invariant 
description with a true dynamics of gauge invariant observables which directly correspond to the usual 
metric and matter fields. 

Instead of Einstein's equations, we now have the Hamiltonian equations for those fundamental fields. 
The interesting question is then how close these equations come to the Einstein equations. In [1] we showed 
that by identifying suitable fields, one can almost exactly match our equations of motion with Einstein's. 
However, there are corrections which are due to the presence of the dust. More precisely, the dust itself no 
longer plays a role at the level of the physical observables, however it implies the existence of a conserved 
energy momentum current density with appropriate conservation laws which substitute the usual initial value 
constraints of General Relativity. This is potentially dangerous since, due to the absence of the constraints, 
our theory has four more physical configuration degrees of freedom than it would have without dust. Hence 
our theory should predict that these modes are either not excited or have decayed in order to agree with 
standard General Relativity (two physical degrees of freedom) plus the observationally confirmed matter 
of the standard model. Of course, it would be intriguing if one would find those additional modes in an 
experiment which then would provide indirect experimental evidence for the existence of the dust. 

How much these corrections and additional degrees of freedom affect the results on the solutions of 
Einstein's equations is difficult to decide in full generality and analytically, therefore we will resort as usual to 
perturbation theory. In [Ij we applied linear Hamiltonian perturbation theory to our system of Hamiltonian 
equations of motion. That is, one takes a general exact solution and then perturbs the equations of motion 
to linear order in the perturbations. Notice that in contrast to usual perturbation theory our perturbations 
are manifestly gauge invariant, they are full (diffeomorphism invariant) observables. 

In this paper we specialise the formulae found in [T] to an FRW background which is an exact solution 
to our equations of motion. This allows us to study linear cosmological perturbation theory. Needless 
to say, this is a topic of fundamental importance in modern cosmology which goes back all the way to 
the pioneering work by Lifshitz [6]. Early on, it was realized that the problem of gauge freedom poses 
serious challenges for the use and interpretation of perturbation theory in general relativity, in general, and 
cosmology, in particular. A significant breakthrough was the construction of quantities that are gauge- 
invariant to linear order [7]. Based on these, a satisfactory framework for linear perturbation theory for 
cosmology was developed and then applied by many authors. Authoritative review articles can be found in 
[U [9]. For the purpose of this paper, we will call the standard approach as presented in these references 
"standard cosmological perturbation theory" and abbreviate it as SCPT. In that context, we will always use 
the notation of [9]. While there seems to be little room left for significant improvements at the linear level, 
developing perturbation theory to higher orders is far from fully established, due to the challenge of finding 
gauge-invariant quantities. Some recent approaches can be found in [10 \ [T2l 1131 H^j. That the question 
of higher order perturbations is far from academic can be seen, for instance, from the growing interest in 
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the issue of non-Gaussianity of cosmological perturbations |TBj. This is where we expect our approach to 
lead to major progress, due to its method of implementing gauge-invariance. Before following this line of 
investigation, however, we should see the results that our framework delivers at linear order and compare 
them with SCPT. 
We find that 

1. We can match our gauge invariant geometry and matter observables with the linearly gauge invariant 
ones of SCPT. 

2. We can cast our equations of motion into SCPT form plus corrections which are proportional to the 
conserved energy momentum tensor current 

3. The corrections decay in the late universe at least as 1/a where a is the physical scale factor and are thus 
negligible whatever their value is. This is crucial because otherwise our theory would predict, even without 
matter, the coexistence of vector and scalar modes apart from the gravitational waves. 

We conclude that, at least in cosmological settings, our theory is in agreement with the usual results. As a 
special case this also extends to the Minkowski background (standard model plus gravity, in vacuum gravi- 
tational waves). Hence our theory has passed a first important consistency check. It should be mentioned 
here that there is an alternative approach to cosmological perturbation theory, based on the use of covariant 
quantities \16\ [T7] . We already compared it with our framework in [IJ , so here we will limit ourselves to a 
detailed comparison between our results and SCPT. 

The plan of the current paper is as follows: 

In section two we specialise the gauge-invariant linear perturbation theory of pj to the FRW background. 
In appendix B, as a calculational check, we also follow a shorter route which is to first linearise the first time 
derivative equations and then to derive the second time derivative equations for the perturbations, which, 
as we know from appendix D of pT], must give the same result. 

In section three we review the SCPT framework from a conceptual point of view. We then define a map 
between the linearly gauge invariant SCPT variables for gravitational perturbations and our fully gauge 
invariant variables. We show that under this map our equations of motion can be cast into SCPT form plus 
corrections which, as already stated, decay. 

In section four we repeat this analysis for the case where dust is included in SCPT. We then show that 
the corresponding SCPT perturbation variables precisely match our gauge invariant perturbation variables 
when we expand them up to first order in the non gauge invariant perturbations. This demonstrates that 
the identification performed in section three is correct. 

In section five we summarise and conclude. 

In appendix A we give a concise review of SCPT from a technical point of view. 

In appendix B we compute the equations derived in section two by an alternative and faster route, thus 
providing a consistency check. As a side result, we show that the ordinary FRW background is a solution 
also to our theory (with a modified energy density) and thus provides a viable background. 

For the benefit of the reader, we list our notation once more on the next page. 
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Notation 



As a rule of thumb, gauge non invariant quantities are denoted by lower case letters, gauge invariant 
quantities by capital letters. The only exceptions from this rule are the dust fields T, , p, Wj, their conju- 
gate momenta P, Pj , /, P and their associated primary constraints Zj , Z, Z^ which however disappear in the 
final picture. Partially gauge invariant quantities (with respect to spatial diffeomorphisms) carry a tilde. 
Background quantities carry a bar. Our signature convention is that of relativists, that is, mostly plus. 
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symbol 



meaning 



Gn Newton constant 

K = IGttGn gravitational coupling constant 

A scalar coupling constant 

A cosmological constant 

M spacetime manifold 

X spatial manifold 

T dust time manifold 

S dust space manifold 

v, p, .. = 0, .., 3 tensor indices on M 

a, b, c, .. = 1, 2, 3 tensor indices on X 

k, .. = 1, 2, 3 tensor indices on S 

Xf^ coordinates on M 

x"" coordinates on X 
coordinates on S 

t foliation parameter 

r dust time coordinate 

Yf' one parameter family of embeddings S ^ M 

St = Yt{S) leaves of the foliation 

g^y metric on M 

Qah (pullback) metric on X 

(jij (pullback) metric on S 

Qij Dirac observable associated to qab 

p"^ momentum conjugate to Qab 

p^^ momentum conjugate to q^j 

P^^ momentum conjugate to Qij 

C scalar field on M 

^ scalar field on X 

^ pullback scalar field on S 

H Dirac observable associated to ^ 

TT momentum conjugate to ^ 

TT momentum conjugate to ^ 

n momentum conjugate to S 

V potential of C, i-, S 

T dust time field on X 

T dust time field on S 

dust space fields on X 

p dust energy density on M, X 

Wj dust Lagrange multiplier field on M, X 
U = ~dT + WjdS^ dust deformation covector field on M 

J = det{dS/dx) dust field spatial density onX 

P momentum conjugate to T 

P momentum conjugate to T 

Pj momentum conjugate to 

I momentum conjugate to p 

P momentum conjugate to Wj 

Zj, Z, Z^ dust primary constraints on X 

p,^ , p, Pj dust primary constraint Lagrange multipliers on X 
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difFeomorphism of X 

unit normal of spacelike hypersurface on M 

n coordinate lapse function on X 

n" coordinate shift function on X 

p momentum conjugate to n 

Pa momentum conjugate to n" 

z, Za primary constraint for lapse, shift 

ly, v"" lapse and shift primary constraint Lagrange multipliers 

0, V', B, E SCPT scalars on X, S 

Sa,Fa SCPT transversal vectors on X 

Sj, Fj SCPT transversal vectors on S 

hab SCPT transverse tracefree tensor on X 

hjk SCPT transverse tracefree tensor on S 

linear gauge invariant completions of ^, ^ 

Va linear gauge invariant completions of Fa 

Vj linear gauge invariant completions of Fj 

c*°* total spatial diffeomorphism constraint on X 

c*°* = SjC^a^ total spatial difFeomorphism constraint on X 

c*°* total Hamiltonian constraint on X 

Ca non - dust contribution to spatial diffeomorphism constraint on X 

Cj = SjCa non - dust contribution to spatial diffeomorphism constraint on X 

Cj non - dust contribution to spatial diffeomorphism constraint on S 

Cj 7^ Cj momentum density: Dirac observable associated to Cj 

c non - dust contribution to Hamiltonian constraint on X 

c non - dust contribution to Hamiltonian constraint on S 

C ^ c Dirac observable associated to c 

h energy density on X 

h energy density on S 

H = h energy density: Dirac observable associated to h 

hj = c*°* — Pj auxiliary density on X 

€ numerical energy density on S 

Cj numerical momentum density on S 

H = d?a H physical Hamiltonian, energy 

L Lagrange density associated to H 

^ ~ fs ^'^^ ^ physical Lagrangian 

Vjk velocity associated to Qjk 

T velocity associated to S 

N = C/H dynamical lapse function on S 

Nj = —Cj/H dynamical shift function on S 

= Q^^N]^ dynamical shift function on S 

V ^ g^u compatible covariant differential on M 

Da Qab compatible covariant differential on X 

Dj Qjk compatible covariant differential on S 

Dj Q jk compatible covariant differential on S 

Qjl^ background spatial metric 

P^'^ background momentum conjugate to Qjf. 

S background scalar field 

n background momentum conjugate to S 

_ • 2 _ 

p = jv[E + v{'E)] background scalar energy density 
_ • 2 _ 

p = ^[S — ■u(S)] background scalar pressure 



Gjkmn = Qj{mQn)k " \QjkQmn physical DeWitt bimetric 

[gf-ip-fcmn ^ QjirnQn)k _ QjkQmn jn^erse physical DeWitt bimetric 

Gjkmn = Sj(m^n)k — \^jk^mn flat background DeWitt bimetric 

p-^^^3kmn ^ §j{mgn)k _ gjkgmn inverse flat background DeWitt bimetric 
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2 Specialisation of Linear Perturbation Theory to FRW Background 



In this section we want to specialise the perturbed equations that were derived in our companion paper 
[1] for an arbitrary background to the case of an FRW spacetime. The perturbed matter scalar field S3, 
as well as the perturbed three metric SQjk are manifestly gauge invariant. This gauge invariance is not 
restricted to linear perturbations only, but extends to arbitrarily high orders. The underlying reason is the 
fact that we applied perturbation theory to an already manifestly gauge invariant physical system, since the 
quantities Qjk and H are gauge invariant by construction. Hence, any function of 6E, 6Qjk is still a gauge 
invariant quantity. For this reason higher powers of the perturbation that occur in higher order perturbation 
theory will in our framework not destroy gauge invariance. Therefore, in principle, it is possible to analyse 
cosmological perturbation theory also in arbitrarily high orders. However, in this section we will focus on 
cosmological perturbation theory up to linear order, the discussion on higher order equations will be left for 
future work. 



Due to the homogeneity and isotropy associated with FRW, the final second order equations of motion, 
shown in equation (6.13) in [1] for the linear perturbation of the matter scalar field 6'E, and in equation 
(6.31) for the perturbed three metric dQjk, simplify drastically. A similar simplification occurs for the co- 
efficient functions that were introduced in section 6.3 in [1] in order to write those second order equations 
in a less complicated form. Before actually discussing the specialisation to FRW for the perturbed equa- 
tions of motions, we will first consider the unperturbed second order time derivative equations of motion 
for S and Qjk, derived in [1] in equation (4.24) and (4.25), respectively, and check that they yield a gauge 
invariant version of the ordinary FRW equation of standard cosmology. An alternative derivation of the 
gauge invariant FRW equation can be found in appendix [B] in section [B.H where already the first order time 
derivative Hamiltonian equations are specialised to FRW before the second order time derivative equations 
are calculated. 



2.1 Gauge Invariant FRW Equations 

In the remaining part of this paper we will indicate quantities which are specialised to FRW spacetimes 
with a bar. Due to the homogeneity of an FRW spacetime, all spatial derivatives vanish. For this reason the 
matter as well as the geometry part of the (gauge invariant) diffeomorphism constraint, denoted by Cj 
and C^^° , vanish as can be seen from equation (j2.1|) : 



C(r,cT) 



[-2Q^,{DkP'') + nE,,]{T,a) 
1 r 



= \QjmQkn 



VdetQ 



-QjkQ ^ 



pjkpmn 



detQ R^^\Q]+2A\/detQ 



r,cr 



1 r 
+ 2A 



H^ 



VdetQ 



+ Jdet Q [g''' (DjE) (DkE) + v{E)] (r, a) 



C'' {T,a)+C (r,cj) 



(2.1) 



where Dj is the covariant differential compatible with Qjj^ and we used that the Christoffel symbols with 
spatial indices vanish in the case of FRW. The results in pQ showed that one of the effects when dust is used 
as a clock is that we obtain a phase space dependent and thus dynamical lapse function and shift vector, 
which we denoted by N and Nj, respectively. Their explicit form is given by 



H{a) 



and 



H{a) 



(2.2) 
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where 



H{a) = ^C{T,aY - Q'\T,a) C,{a) Ck{a) (2.3) 



is the physical Hamiltonian density. This is the density of the physical Hamiltonian H that is generating 
the (dust) time evolution for all observables, in particular for H and Qj^- Knowing now that Cj vanishes 
for FRW, we get immediately that the shift vector vanishes, too. Hence in the case of FRW, all terms 
proportional to Nj can be omitted in equations (4.17) and (4.18). Furthermore from equation (j2.3p and the 
definition of N and Nj we get for FRW 



N= = = ^ = — = \ll + Q' N,Nk = l. (2.4) 

We observe that N and N j are not independent quantities and that for the special case of FRW the lapse 
function takes the constant value +1. Keeping this in mind and taking into account that all terms containing 
spatial derivatives vanish, the equations (4.17) and (4.18) for H and Qjk in our companion paper [1] reduce 
to 

^ ^ (v^d^y ^ 1 

7d^ ^ 2"^^^^ 



= - ^"^^ + + ( - -7f^ C + 2A+ ^v{E) ] . (2.5) 



Using that Q^^ = A^5jk where A is the physical scale factor, that is the gauge invariant extension of the 
usual non gauge invariant scale factor a used in standard cosmology, we obtain for H 

'i\= 1 



This agrees (formally) with the usual FRW equation for a scalar field. The agreement is formal in the sense 
that here we have an evolution equation for a gauge- invariant and thus physical scalar field H. Also the 
equation contains the physical scale factor A, not a which usually occurs in this equation. Furthermore the 
dot refers to a derivative with respect to dust time r in the equation above. For Q^^ we need the explicit 
form of the matter, as well as the geometry part of the Hamiltonian constraint, whose sum is denoted by 
C. Specialised to FRW it is given by 



Inserting this back into equation (|2.5p and performing the derivatives yields 

i)=-2(a) +2'^- 4A2(= <^-**> 

As shown in our companion paper yp, the Hamiltonian density H in equation (j2.3p is a constant of motion. 
In the case of FRW, H reduces to 



H = Jc'^ - Q'^NjNk = \l& = C (2.9) 



due to the vanishing of the shift vector. Denoting the energy density by e, we have the conservation law 

H = C = A^ 
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Note that C is not a constraint here, since only the total sum C*°* = C^'^" + (j"^^^^^"^ _|_ with the dust 



included, is constrained to vanish. The conservation law above can be equivalently written as 



A 



en 



2^3 +.(-)). 



^2 



Replacing [A/A)"^ in equation (j2.8p by the expression in the equation above, we end up with 



^-4 



2AV 



.-(- 



e 

A3 



Introducing the matter density and pressure as well as the dust density 

^^matter 



—matter 



-^ust 



—matter 



c 



A3 

1 /=2 

2A 



-:=7matter -p^matter 



A3 



c_ 

I3 



e 
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we can rewrite equation (j2.12p as 



<5 



—matter , —dust , Q7=rmatter \ 
P 'T~ P -r op I . 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



This again coincides formally with the usual FRW equation for the gravitational part if the system "grav- 
ity + scalar field + dust" is considered and the dust is assumed to be pressureless as is the case in our 
framework. The minus sign in p^^^^ in equation (j2.13p reflects the phantom nature of the dust which was 
discussed in detail in section 2.3.4 of our companion paper [1]. 



Summarising, when (phantom) dust is employed in order to deparametrise the constraints of General Relativ- 
ity, the resulting physical Hamiltonian generates equations of motion that formally agree with the standard 
FRW equations. For this reason, the dust clock seems to be the clock closest to the framework of standard 
cosmology where the Hamiltonian constraint is taken as a true Hamiltonian. In the next section we will 
move on to the discussion of the perturbed equations of motion for 5E and 5Qjk- 



2.2 Specialisation of the Perturbed Equations of Motion to an FRW Background 

Let us now specialise the perturbed equations for 6E and 6Qjk, derived in our companion paper [1], to 
the case of an FRW background. As discussed in the last section, in this case the lapse function takes the 
constant value = 1 and the shift vector Nj vanishes. Taking into account that all spatial derivatives of 
background quantities (these are indicated with a bar) vanish, the coefficient functions for 5E that can be 
found in equation (6.33), (6.34) and (6.35) of [l] simplify to 



[Cs] 



(Vdetgy d 



+ 



VdetQ 
2 dr^^ ^ 



9r2 +^(yl)9r 
1 ^A\^ 



dr 



A'^KA 



1 1 ^..^ d 



2^2 



dr 



dx^ 



1 ^pk 9 

dx^ 



(2.15) 



Using these coefficients, we obtain the following perturbed equation 



r d 



Idr^ KAJdr dx^ dx'' 



1 



d d 



+ 



A^\A 



1 1 ^..i. d 



2^2 



dr 



6Qjk + 



2A^ . 



6N,. 
(2.16) 
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When performing the derivatives and rearranging the terms, the equation above can be rewritten as 



+ ^5^'6E,, - -v"{E)6E + -t^-)E6Q,, - -j^E6Q,, + -^E6N,,. (2.17) 

This equation agrees with equation ()B.22p in appendix [Bl where an alternative route is taken to obtain the 
second order equation of motion. There we first perturb the Hamiltonian first order equation and derived 
then the second order perturbed equations. For hnear perturbation theory it is proven in appendix D of 
[l] that these two ways of deriving the second order equation of motion are equivalent. Compared to the 
perturbed equations in SCPT we obtain an additional term, namely the last term on the right-hand side 
of equation (I2.17p . This term arises from the interaction of the clock, namely the dust, with the physical 
system. However, in the next section where we compare our results in detail with the results from SCPT, 
we will show that the physical predictions of both frameworks are in agreement. 

Let us now consider the coefficients for 6Qjk in equation (6.37)-(6.41) of [Ij and their specialisation to FRW. 
Taking again into account that = 1, that all terms proportional to Nj can be omitted and that all spatial 
derivatives of background quantities are zero, we obtain 

[^q] = TT^ + , = C-{2A+ :rrv{^)) - DmDnQ + 



dr^ 2Vdetg ' 2A ^ y^detg dr 



I K d Ik 



A^2X dr A'^AX 



5. 



'jk 



im _ 7) d—mn (\/det Q) (-^ B —mn d—mn\ 1-^ -ir ^ 

\ Q j^g^ ^ 1 ^ ^ ^ ^ \ ^ ^ Q 

[GQ]jk ~ ~'2^i^'Q^^ ~^ ^ QrjQsk + ^QjkQrsQtuQ [G ] — -^QjkQrs[G ] ^ 

1, 



+DjDkQ""' + -[G-Y'''''DrD, 



5^ + D.DkQ^'' + i[G-i]™Z),Z),. (2.18) 



For the derivation of [Cq]^" we used that i?™-" vanishes in the case of FRW and we also employed several 

times the identity Qmn[G~^]^^"^^ = —2Q^^. Recall that the coefficient [BQ]jk is associated with the pertur- 
bation of the last term on the right-hand side of the equation of motion for 6Qjk in equation (4.25) in [T]. 
We mentioned already that this term is the only deviation from the standard Einstein equations that use 
the Hamiltonian constraint as a true Hamiltonian. However, since this term is quadratic in the shift vector 
Nj, it vanishes for both the background equations and the linear perturbations, when specialised to FRW. 
Nevertheless, we will see below that we still get small deviations from the standard treatment due to the 
shift vector being dynamical. Note that since N = y^l + Q^^NjNj. we get 

SN = - — ^^SQjk + N^6Nj = 0. (2.19) 
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Hence, the perturbation of the lapse function also vanishes in the context of an FRW background. Finally, 
we use the coefficients in equation ()2.18p and obtain the following equation: 



9r2 



— -3 - -A+— fH -V 



A) 



4A 



)) -Z)^Z)„Q"" + 3(-)— ],5g,fc 



d 



3 Q^k 



dxi) 



<5H + 



d 



5T^ + 5l 



dx'' ^ dxi J dr 



A\d_ 
AJd^ 

d _ ^ 
5r 

d \ d 



(2.20) 



6Ni 



-limnrs 



DrD, 



6Q, 



Applying the derivatives to the perturbation terms, rearranging them and using equation (j2.8p . we end up 
with 



Qjk 



A 



A 



A 



2A 



1 



MH)+2(4)5iV, 



A) 



(2.21) 



+DmDnQ 5Qjk + D.DkQ + -[G-T'^^'DrDs 5Qmn " 2DnD^kQ SQj) 



In the equation above we used that the perturbed shift vector 6Nj is a constant of motion which was 
shown in our companion paper p[] in appendix E. Therefore we have 5Nj=0 and thus can omit this term 
in equation 12.201 This equation agrees with the one for 5Qjk in equation (|B.22p if we consider that the 
perturbed Riemann tensor can be expressed as 



6R 



■jk 



and the perturbed Ricci scalar reduces for FRW to 



(2.22) 



(2.23) 



Likewise, in the case of the matter perturbation we obtain an additional term containing the variation of the 
shift vector which does not occur in the SCPT analysis. This finishes our discussion of the specialisation of 
the perturbed equation of motion to an FRW background. In the next section we will compare our results 
in detail with the one obtained in SCPT. 



3 Comparison with SCPT 

The aim of this section is to show that our invariant Hamiltonian perturbation theory reproduces the SCPT 
results up to small deviations. By this we mean the following: 

1. SCPT Lagrangian Approach: Linearly gauge invariant, linear perturbations 

SCPT uses the Lagrangian formulation. It takes as matter a scalar field C coupled to the metric g^u. 
These eleven fields are of course not invariant under (infinitesimal) spacetime diffeomorphisms of the 
manifold M which are considered as gauge transformations. However, one constructs seven functions 
which are certain linear combinations built from the metric perturbations 

6goo = 2a^(j), 6goa = a^iSa + -B,a), Sgab = a^[2Safe + 2ip6ab + 2F(„^b) + hab] (3.1) 

and from the scalar field perturbation 6C. Here a is the scale factor of the spatially flat FRW back- 
ground, the vector fields Fa,Sa are transversal with respect to that flat Euclidean spatial metric and 
the tensor h^b is transversal and tracefree. The seven invariants are explicitly given by 

^:=^-n{B-E')-{B-E'y, ^ :=iP + n{B-E'), Va:=Sa-Fa, hab, Z := 6C + (' (B - E') (3.2) 
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where TC = a' /a, a prime denotes derivation with respect to conformal time dtj = dt/a and Q is the 
background scalar field. 

One then expresses the ten perturbed Einstein equations directly in terms of these seven invariants. 
Four of these equations, namely the temporal - temporal and the temporal - spatial equations do 
not contain second time derivatives of four of the seven fields, they are constraints. They imply that 
four of the seven invariants can be expressed in terms of the other three. The system has only three 
independent, linearly gauge invariant degrees of freedom. Specifically, the transversal part of the 
temporal - spatial part of those equations imposes 14 = 0, while its longitudinal part together with 
the temporal - temporal equation allows us to express, for instance, ^l', Z in terms of ^ . Namely, since 
the cosmological constant term and the energy momentum tensor have no contribution of the form 
f^ab it follows immediatelj0 from the third equation in ()A.54p that <I> = From the longitudinal part 
of the temporal - spatial equations we then immediately find that Z + 4A(^' + 7i^)/K is a spatial 
constant which then must vanish due to the boundary conditions. 

Thus the system of a scalar field coupled to the geometry has only the three degrees of freedom 
hab, ^ as independent, linearly invariant, propagating degrees of freedom. Without Z we even have 
^ = ^ = 0. 

2. Hamiltonian Approach: Manifestly invariant, linear perturbations 

Our approach is completely different: first of all we work in the Hamiltonian framework. Secondly, 
we add additional four dust scalar fields. These are altogether fifteen configuration degrees of freedom 
to begin with and they are accompanied by the corresponding canonical momenta. Thus, from the 
very outset the number of degrees of freedom is very different compared to the SCPT framework. The 
components goo, goa of the metric have the same information content as lapse and shift functions n, Ua 
respectively. Their conjugate momenta are constrained to vanish by the primary constraints, hence 
n, Ua are pure gauge and can be considered as Lagrange multipliers in the Hamiltonian formalism. 
This is in agreement with the Lagrangian formalism where dgoo = 2a'^(j), 6goa = a^{B^a + Sa), which 
are also pure gauge and are used to construct the gauge- invariant functions ,Va, hab- Thus, the 
number of physical configuration degrees of freedom is reduced by four. In addition there are the 
secondary Hamiltonian and spatial diffeomorphism constraints. These reduce the number of physical 
configuration further by four. Thus, we end up with 15 — 8 = 7 physical degrees of freedom. Without 
dust we also would have 3 physical degrees of freedom only, which shows that with equal matter 
content the Hamiltonian and Lagrangian framework are equivalent. 

Indeed [19] the temporal - temporal and temporal - spatial components of the Einstein equations are 
proportional to the Hamiltonian and spatial diffeomorphism constraints, respectively, and the gauge 
transformation equations (in second order time derivative form) generated by the constraints are 
equivalent to the remaining spatial - spatial Einstein equations. The gauge transformations generated 
by the secondary constraints of the variables different from lapse and shift and their conjugate momenta 
(the secondary constraints, when written in terms of the momenta conjugate to the spatial metric and 
the scalar field, do not depend on them) are equivalent to spacetime diffeomorphisms, on shell [l9j. In 
fact, this statement can be extended also to the lapse and shift fields [2QJ as we will recall in the next 
section. 

Hence the difference between our framework and the usual one is not in the usage of the Hamiltonian 
versus the Lagrangian formalism but rather in the addition of the four dust fields. Now what one 
could do is to consider the dust - scalar - geometry system, perform linear perturbations in non gauge- 
invariant quantities and then construct linearly gauge-invariant quantities from those. This we will do 
in the next section. However, what we will do here is to construct quantities which are fully gauge- 

^We assume here that the perturbations vanish sufficiently fast at spatial infinity which is consistent with making asymp- 
totically FRW boundary conditions. Then the Laplacian A has no zero modes and is therefore invertible. An equation of the 
form f^ab + gSab ~ theu implies A/ + 3g = and A(A/ + g) — that is f = g — 0. 
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invariant to all orders and then look directly at perturbations of those manifest invariants. Hence, 
gauge-invariance is treated here non-perturbatively, and our perturbations themselves are manifestly 
gauge- invariant. This is achieved by means of deparametrisation in terms of the dust degrees of 
freedom, which in one stroke solves the secondary constraints and spells out the true i.e. gauge- 
invariant degrees of freedom. Namely, the dust momenta can be expressed in terms of those invariants, 
and the dust fields themselves together with the remaining fields form non linear expressions which 
are fully gauge- invariant. This leaves us, before performing perturbation theory, with seven degrees 
of freedom, namely the invariants H, Qji^ corresponding to the scalar field C (or rather its pull back 
^ to the spatial manifold X) and the spatial - spatial components of the metric together with the 
corresponding canonical momenta. 

In terms of these seven degrees of freedom there are no longer any secondary constraints, they have 
been reduced. However, our equations of motion for the perturbations (5H, 6Qjk of the seven perturbed 
manifest invariants take a form almost completely identical to those equations that we obtain in the 
Lagrangian formalism for the seven perturbed linear invariants ^,^,Va,hab, when proper identifica- 
tions are made. The crucial difference arises from the fact that the latter quantities are still subject 
to constraints while the former are unconstrained. Yet, it turns out that while in our formalism there 
are no constraints any more, the would - be constraints are now constants of the motion. These 
constants of the motion express the influence of the dust on the system. In particular, in the limit of 
vanishing influence (test dust) these conservation equations turn effectively into constraint equations 
and then we get a precise match between the two formalisms with proper identification of variables 
and spatial manifolds X, S, respectively. Thus, for test dust our formalism just produces arbitrarily 
tiny modifications, that is, source terms, of the usual formalism although fundamentally the number 
of physical degrees of freedom was changed. In particular, in the late universe these modifications 
decay as compared to the usual terms. 

It is quite remarkable that the two formalisms with fundamentally different numbers of degrees of freedom can 
be brought to such a close match. Looking at the details, it is hard to imagine that any material reference 
system other than pressure free dust can lead to such a modest modification of the usual formalism. Indeed, 
in the appendix of [1] we compute the modifications that arise when using instead the phantom scalar 
field of [5] for purposes of deparametrisation. It is shown that the modifications, in fact, grow in the late 
universe which is qualitatively drastically different. However, the coupling constant of the phantom and 
the constants of the motion can be tuned in such a way that the moment of time when the modifications 
become significant can be chosen to be arbitrarily late. 

It is worth mentioning that since we can take the limit of 1. vanishing dust influence (i.e. both the 
background energy momentum (e, Cj) and its perturbation {6e,6ej) can be taken to zero), 2. vanishing 
cosmological constant, 3. vanishing physical scalar field and finally 4. vanishing time dependence of the 
background FRW scale factor a = 1, it follows immediately from our formalism that in this limit we recover 
also the two gravitational wave modes on Minkowski space, the other four modes are frozen out by the 
equations of motion. 

The reader may now ask what the advantage of our more complicated formalism is compared to the usual 
one. There are at least three worth mentioning: 

A. Higher order classical general relativistic perturbation theory 

As far as standard cosmological perturbation theory is concerned, the usual formalism cannot be eas- 
ily extended in a gauge invariant way to higher than linear order. The reason is that one has to 
construct the relevant gauge-invariant quantities at each order from scratch again. In our manifestly 
gauge- invariant formalism, by contrast, those quantities are simply given by the relevant higher-order 
perturbation of the metric and matter degrees of freedom. This is a tremendous simplification. We 
expect that this should also make it much easier to address general question of the stability of cosmo- 
logical perturbation theory [23] . 
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B. Conceptual Improvement 

General relativity is a generally covariant field theory without true Hamiltonian. This leads to many 
conceptual problems such as the problem of time (gauge invariant functions do not evolve with respect 
to the Hamiltonian constraint (s)). The reason that in perturbation theory one still has non trivial 
evolution of the linear invariants is only due to the fact that one has chosen a background spacetime 
which provides a preferred notion of time and what one computes is evolution with respect to that 
background. However, fundamentally general relativity is a background independent field theory and 
since when going to higher order both classical perturbation theory (due to reasons of gauge invariance) 
and quantum perturbation theory (due to reasons of non renormalisability) fail, one is forced to adopt 
a non perturbative approach so that the problem of time reappears. The relational formalism solves 
the problem of time by determining a physical Hamiltonian (which is not constrained to vanish) from a 
chosen material reference system which evolves the manifest invariants. This brings the analysis back 
to the conceptually safe realm of a dynamical Hamiltonian system, albeit a technically complicated 
one. 

C. Quantum Theory 

At least for the pressure free dust chosen here, the Poisson algebra of the linear invariants is as simple 
as for their non invariant analogs. Therefore it is possible to adopt a reduced phase space quantisation 
approach [2T] and to find Hilbert space representations of that algebra. This has the advantage that, 
for instance, in the framework of Loop Quantum Gravity the kinematical Hilbert space used there 
becomes now a physical Hilbert space. The constraints have completely disappeared from the screen. 
What remains is to quantise the complicated physical Hamiltonian [22] and to analyse its spectrum. 
This may be technically complicated but, conceptually, it is crystal clear. In particular, the fact that 
one has a physical Hamiltonian at one's disposal which by construction does not depend explicitly 
on some background time may improve the vacuum problem that one encounters in QFT on curved 
(time dependent) backgrounds. The problem with explicitly time dependent Hamiltonians is that the 
notion of vacuum (ground state) is time dependent which may lead to (infinite) particle production 
and, in particular, means that there is no unique ground state. The physical Hamiltonian, in principle, 
selects a preferred class of states as ground states, namely its zero eigenvalue states. It may be that 
zero or the minimum of the spectrum is not in the discrete (more precisely, pure point) spectrum or 
that the zero eigenvalue is vastly degenerate, however, at least conceptually, this appears to be an 
improvement. 

In what follows we now turn to the proof of the above claims. 
3.1 Conservation Equations 

As shown in the appendix of , for any fully conserved quantity -F of a Hamiltonian sytem with Hamiltonian 
H, when we expand both the equations of motion and F to order n then F is still a constant of motion up 
to terms of order n + 1. In particular, for n = 1 the coefficients of the correction vanish. This means that 
we can derive conservation equations for manifestly gauge-invariant perturbation theory by expanding the 
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four times infinitely many constants of motion Cj{a), H{a) where 



H{a) 
C 



1. 



—DkP^ + -liD,~ 
^C7(a,r)2-gjfe(a,r)C,((T)Cfc(a) 

C'geo + C\ 



matter 



geo 



1 



matter 



2A 



+ x/dit(Q)(Q^''SjH,fc + v{E)) 



(3.3) 



Let ej(o') = — Cj(cj), e{a) = H{a) be the zeroth order values of the constants of motion which are obtained 
explicitly by inserting the FRW solution of appendix IB. II into the expressions ()3.3p . We now expand 
5ej := ej — ej, 6e := e — e to first order, which then are still constants of the linearised equation of motion. 
We find after a short calculation (remember that DkPj = t^fc-Pf ~ ^kj-^ t^ since P^^ is a tensor density) and 
using the perturbed gauge- invariant variables defined in appendix IB. 21 



6e = = 



,k 



--{A'SPf -2A{6Q 



K(5Cgeo 

6Rjk 

A^Cjnatter 



C5C - CkSCj - CkCj5Q^^ 

i2 



A 

5C 



—6^^6Qjk - 2AA6jkSP^'' + Ad^^SRjk + kAP^SQjk 



1 



2A^ _ 
1 r IT 



2SQi{j,k)i - ^^Qjk - ^Qu,jk 
1 



2^6U - -AW'SQ.k + A'v'{r)5E 



(3.4) 



Here we used that that the quantities L,I introduced in appendix lB.il are explicitly given hy L = and 
/ = —2 A. We will discuss the momentum and energy conservation equations separately in what follows. 



3.1.1 Momentum conservation equation 

We substitute the perturbed configuration variables for the perturbed momenta by using the equations of 
motion, see equation (jB.2ip . Notice that Nj = —Cj/H is also a constant of motion and that 5Nj = 5ej/e. 
Then we find from the first relation in (j3.4p 

kA 

6ej + — [A5ej - 5ek,kj] 
= - [A{6Qjk^k - SQkk,j) - 2A{5Qjk,k - SQkkj)] - T-ntJHj. (3.5) 

K A 

In order to connect (j3.5p with the SCPT equations we parametrise our metric perturbations in the following 
way: First of all, SCPT employs conformal time = t] for which ^oo = — a^- On the other hand, we always 
have Qrr ■■= -N^ + Q^^NjNk = -1, which means that we automatically work with cosmological time, that 
is, the eigentime of the dust. We use the letter A instead of a for the scale factor in what follows in order to 
emphasise that A is an observable, rather than a gauge-dependent function. From the relation dx^ = dr/A 
we find Qrr = 500/^^) drj = 5oj/^ = ^j-, 9jk = Qjk- Now in analogy to SCPT we introduce the following 
ten functions 

500 = (-1 + 24>)A^, goj = (Bj + 5,)^^ gjk = A^{2i^5jk + 2E^,k + 2F(,- fc) + hjk), (3.6) 
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where the tensors Sj,Fj,hjk are transversal with respect to the Euchdean metric 5jk and hj^ is trace free. 
We can also express ip, E, Fj , hjk as 



4^2 
1 



k 



Fj = -^A-^6Qjk,k- A'^SQkijki) 

hjk = 6Qjk-2{'^6jk + Ejk + F^j^k)) ■ (3.7) 

From g-TT = — 1 we immediately see that automatically 

<A = 0, (3.8) 

that is, we are forced to work in partly synchronou^ "gauge". We write this expression in inverted commas 
because there is no gauge involved here, all the quantities are manifestly gauge invariant. We use this term 
only to make the analogy clear. In particular, since Nj = —Cj/H can be expressed in terms of Qjk, ^ 
and (via the equations of motion) their velocities, in our formalism the functions B, Sj are not independent 
variables from the outset. 

Let us decompose (j3.5p into longitudinal 5e^j = A^^dck^kj and transversal dej- = Cj — e'-' parts, respectively 



. ^ , 1 



64 + '^^ASej 



A {5Qjk,k - A~^6Qki,kij^ - 2A [5Qjk,k - A~^5Qki,kij) 



Then we find 



6ei 7 = — 



A [sQjk,jk - ASQkk) - 2A {5Qjk,jk - AdQkk)] - ^HA^S. (3.9) 



Next we invoke dMI), (I2Z3) remembering that 5Nj = A{Sj + Bj) = -5ej/Z and that (!) = d{.)/dT = ^(.)'. 
Then the transversal part of p.9p becomes simply 

AV,- = -K-^ (3.10) 

where we have introduced the variable Vj = Sj — -Fj which is the analog to the variable Va in the SCPT 
framework, see appendix [Al We also used /' = af several times and, of course, L = and I = -2A. 
Equation (j3.10p should be comparecH with the vector contribution of the perturbation of the temporal - 



■^In the usual not manifestly gauge invariant formalism, synchronous gauge means (j) = B = while longitudinal gauge means 
B = E = 0. 

■^In order to avoid confusion of the reader we point out that the individual contributions of geometry and matter to the 
Hamiltonian and spatial diffeomorphism constraints (without dust) c and Ca respectively are not exactly the temporal - temporal 
and temporal - spatial components of the Einstein and energy momentum tensor respectively, rather we have the identities 



^^gr-av _ _,^(jgt(q,)x'^n''G^y, Kc®''''" = —2y^ det{q)n'' n'^ G fil, (3-11) 

where n'' is the unit future normal to a foliation X. Noticing that n'' = (X*^ — n"X'^)/n and that Sua = 5{qabn'') = a{B^a + Sa) 7^ 
one finds 

KSca = -a^[SGoa + {B,a + Sa){2H' fic = -2a[SGoo + SH^S^P + AE)] (3.12) 
The additional terms cancel when we replace G,i^ by _E^,^ = G^^ + Ag^i, — ^T^i, since Efj,i, — 0. 
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spatial part of the Einstein equation, written in hnearly gauge-invariant form, which we derived in the 2nd 
equation of ()A.42p . the second relation of ()A.27p and the first relation of ()A.50p . That is 

= 5GQa + ^^gOa - -^^^Oa 

-2(M/' + n<^)^a - \AVa - {2n' + n^)V^ + AaV, - (^(C'^,a + a^ma)) 

-2(M/' + H^),a - ^Ak) - (|^C'^,a) (3.13) 

where we used the spatial - spatial part of the background FRW equations ()A.14p . (jA.lSP in the second 
step, that is 

2n' + 71^ = a^ [k-'^p] (3.14) 



2 

_ ^2 _ _ 

and where p = -^{r. — v{'S)) (this holds with or without dust). The transversal piece of (|3.14p is simply 



- ^AV; = 0. (3.15) 

If we compare this to ()3.10p and if we identify the SCPT variable Va with Vj then we see that, instead of a 
Laplace equation without source, we have a non-vanishing source. However, recall that 5ej' is a constant of 
motion. Hence, the source term decays as and thus plays no role anymore in the late universe. This 

happens the sooner the smaller 5ej-. 

Next we turn to the longitudinal part of (|3.9p which can be written as 

= -I {-^M + r'^-.) (3-16) 

where we used H = A^r, = A^r! . A direct comparison between (13.16P and SCPT should be in terms of 
variables which are direct analogs of the linearly invariant SCPT variables <I>, Z. These are given by 
^1 = ip-\- Tif, ^ = (p — Tif — /', Z = + f where we of course identify C with our H and where f := B — E' . 
Hence we should rewrite ()3.16p in terms of these variables, recalling that in our formalism </> = so that our 
$ actually reads $ = -Hf - f. This leads to 



4A 



In order to arrive at this form of (|3.16p . we simply have rewritten ip,'^ in terms of ^, Z and made use of 
the gauge- invariant background FRW equations with dust (]B.10p . (jB.lip . They are an exact solution to 
our invariant Hamiltonian equations and allowed us to write k(H )^/(4A) = — (W — Tl'^ — eK;)/(4yl) so that 
certain terms really combine to We now compare this to the longitudinal part of the linearly invariant 
projection of the temporal - spatial part of the Einstein equations (j3.13p in the SCPT framework which can 
be rewritten as 



K 



"if' + n^ + —C z =0. (3.18) 

4A J ,a 

We see that (I3.18P and (j3.17p agree with each other under the made identification C = H up to the corrections 
on the right hand side of (j3.17p . The first of these terms decays in time as l/A^ since Se^j is a constant 
of motion. The second term decays only as 1/A under the assumption that B — E' remain small during 
the evolution, an assumption that one always makes in perturbation theory. Notice that in the SCPT 
framework the quantities i?, E and also the combination B — E' is not linearly gauge-invariant. However, 
in our framework, all those quantities are gauge invariant to all orders from the outset so that (I3.17P is a 
consistent relation among manifestly gauge- invariant variables. We just have written it in a form as close 
as possible to the equations in the usual formalism. 
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3.1.2 Energy conservation equation 

Finally, we consider the second conservation equation in (j3.4p which after some algebra can be brought into 
the form 



6a 



matter 



^ (2EdE + v'iE)5Ej + (3^ + A^) p 



(3.19) 



_ -.2 _ 

where p = -^{E + t'(H)). Using our FRW background equations (jB.lOp and (|B.lip . we can simplify (|3.19p 
to 



/ — — ■ 



AA 



-A^ + sn'ip') 



5e - eiSil^ + 6E) 

which now needs to be written in terms of the SCPT-like variables 

*:=V + W/, ^ = -nf-f', Z:=E + E'f with f = {B - E'). 



(3.20) 



(3.21) 



We do this by substituting ip^E in (j3.20p by ^, Z, respectively, thereby picking up the following correction 
term just involving / 

-nnfY - ^{2E'{E'fy + A\'{E)E'f) 



A 



2A 



We want to massage this into a form that involves <I>. We have, using A'^p = ^((h')^ + v{E)), 



= {A'pyf + 2{A'p)f'-^{Afyv^ 

= {A'pyf + 2{A'p)f' + ^'^v{E) 
where we used $ = —{Af)' /A. Now we use our background FRW equation in the form 



(3.22) 



(3.23) 



sn^ = A^ 



A + 



K ( p — e 
~A^ 



(3.24) 



from which we can compute (A'^'p)' by taking its derivative. Combining then (j3.23p and (|3.22p we see that 
the ^'H{T~ify term in ()3.22p is cancelled and that the remaining terms are either proportional to or to e. 
The end result is 

^ [6e-e (3V' + AS - 2{B - Ej + n{B - E'))] 



A 
1 

2X 



2E'z' + A^v'{E)z) - - (-A^ + 3n'^') - A^ 



A K, J 



(3.25) 



This should now be compared with the temporal - temporal Einstein equation, written in terms of the 
SCPT variables, derived in the first equation of (IA.42p . the first equation of (1A.43P and (IA.54p . that is 



= 5GoQ + k5~goQ - -5Too 

= 2(-A^ + m^') + 2a^k^ - ^ (2c' Z' + a^v'(C)Z - 2a^v(C)<^ 

4A V 



(3.26) 



Comparing (j3.25p and (j3.26p . we see that the two equations agree, with proper identification of the variables, 
up to the correction term on the left hand side of (I3.25P which is proportional to 6e, I and the perturbations 
B, E, -0, which in our formalism are gauge invariant. Assuming that these terms stay small during the cosmic 
evolution as one always does, the correction term decays as 1/A and thus is negligible in the late universe. 
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3.2 Evolution Equations 

We derive the following perturbed evolution equation in (jB.26p : 

A ■ A A 

^Qjk = -^SQjk + '^^{^N(^j)k) + 2—6Qjk-26Rjk 

+Sjk (^\dRmm - - ■ (3.27) 

First of all, we change the time variable again to dx^ = dr/A and write dQjk = A'^Hj^, which gives 

6Qjk = A[2HHjk + H'^k) 

5Rjk = 2 {'^^Khk)i — ^Hjk — Hiijk) ■ (3.28) 
Inserting iK28\\ into (K27\\ and using 5Nj = A[Bj + Sj] yields 

2nHj^ + Hj^ = 2n [5 (j + 5 (j] - 26Rjk 

+6jk (^IsRm^^ - ^A\E6E - lv'{E)6E)^ . (3.29) 

Equation (j3.29p should be compared to the spatial - spatial part of the Einstein equations written in linearly 
gauge invariant form. Using the third relation in (|A.42p . the third relation in ()A.43p . as well as the spatial 
- spatial components of (|A.45p . (jA.SOp and (|A.54p gives 

6Gab + A~gab - '^Sfab = \ {K, + 2W/i,b - AKb) - {2n' + n^) Kb - (^(',,6) + 2Hy(„,,) 

+ - ^],afe - [H^ --^) + 2 + {2n' + + ^) + 7^(2^ + $)')] 5ab 

+Aa2 {hab + 2^5ab) 

~JX {^P^'^^-b + 2\pa^^Sab + ^(2C' + Z') - v'(C)a^Z)6ab^ 
= \ Kb + '^'Hhab - AKb) - {V(^^,) + 2nv^a,b)) 

+ - ^],afe - [^(^ - ^) + 2 + {2n' + n^)<^ + ni2^ + $)')] Kb 

^ (2C' ($ + Z') - v'iOa^z) 6ab. (3.30) 



4A 

Here we used p.l4p twice in the second step. 

In what follows we decompose (I3.29P into the various irreducible pieces. We will see that, in contrast to 
the conservation equations, maybe not very surprisingly, the evolution equations do not adopt any correc- 
tions. 

3.2.1 Tensor contribution 

In this case Hjk = hjk, 6Nj = and thus 6Rjk = —Ahjk/2, so that the tracefree, transversal contribution 
to (|3.29p reduces to 

- Ahjk + 2nh'jk + h]k = 0. (3.31) 



21 



This should be compared to the corresponding tensor contribution to (j3.3Up which is given by 

^ [Ki, + 2nhab - Ahab) = 0. (3.32) 

Thus we obtain an exact match upon identifying /ij^ on the dust space S with hat on X which is diffeomorphic 
to 5. 

3.2.2 Vector contribution 

In this case Hji^ = 2A'^F(^j^i^-j and B = 0. A short computation reveals that 6Rjk = 0. Hence, dropping the 
term proportional to 6jk in (j3.29p which is a scalar contribution, we find 

4WF^.,) + 2F^^,) = 2H%fc). (3.33) 

Using the conservation equation (SNj)' = (ASj)' = 0, we find the identity 

nSj + S'j=0 + = 0. (3.34) 

Adding twice the second zero in (j3.34p to (j3.33p and dividing the resulting equation by two, we find 

2W(,-fc) + V(;.,) = 0, (3.35) 

where we have again introduced Vj = Sj — F'-. 

On the other hand the vector contribution to (j3.30p is obviously 

^(U)+2W^Kfe)=0' (3.36) 

which is again an exact match when identifying Va and Vj. It is instructive how the conservation law (j3.34p 
found its way into the equations in order to establish this match. 

3.2.3 Scalar contribution 

This time Hj^ = 2{'tp5jk + Ej/.) and Sj = 0. A short calculation reveals that 6Rjk = —{ipjk + ^i^^jk) and 
thus 5Rkk = —4:Aip. Thus (j3.29p becomes 

m{ij'6,k + E[^k) + 2«<5,fe + E';^^) = 2nBjk + 2(^jfe + Ai;6jk) + 5jk (-2AV - " ■ 

(3.37) 

Equation (j3.37p is of the form f^jk + g5jk = 0. Taking the trace and operating with djdk, respectively, we 
learn that ?>g + Af = and A(q + A/) = 0, respectively. Since, due to our boundary conditions, there are 
no zero modes of the Laplaciaro, we conclude f = g = 0. Therefore we can study the {.)5jk and (.)ja: pieces 
in (|3.37p separately, which we note as 

{m^P' + 2< - 2AV^) 6jk = (^-2Aij - ^A\t5E - ^v'{E)5E)^ 6,k 
[AHE' + 2E" - 2nB - 2^ = 0. (3.38) 



*The proof is standard: Assuming that / is an at least twice difFerentiable function which decays at infinity at least as 1/r 
where r is an asymptotic radial variable. It follows that /,-, is square integrable with respect to L2{S,d^a). Now suppose that 
A/ = 0. Then =< f,Af >= — where the boundary term at the sphere at infinity drops out since //,-,■ decays as 

1/r^. It follows that fj = a.e. and since / is in particular continuous, it follows / =const. But that constant must vanish 
since / decays. 
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Consider first the second equation. We note the conservation equation 6Nj = (ABj)' = 0, which can be 
written in the form 

HB^j + B'j = HB + B' = 0. (3.39) 
Here the second relation follows again from the fact that the Laplacian has no zero modes. It follows that 



$ := - [n{B - E') + {B- E')'] = HE' + E" . 



(3.40) 



Now also using the definition ^ = ip + 7i{B — E') and (j3.39p . a short calculation reveals that the second 
equation in (|3.38|) simply becomes 



Next we consider the first equation in (j3.38p . which can be simplified to 



(3.41) 



(3.42) 



We write SE = Z — r! f , where f = B — E' , so that on the right hand side of (j3.42p 6E is replaced by Z plus 
the correction term 



4A 
4A 

K 

4A 
4A 



{[E'fyf + 2iEyf'-A\viE)yf 



A\{E)] f + 2{E'rf' + 2A'nv{E)f 



2A {A^p)'f + 2{E')\f' + nf) - 2-H ( (h')^ - A\{E)) f 
2\ {A^pY f - 2(h')2$ - AXHA^pf 

f 



2A^ ^ 2 



[A^p]' - 2nA^p 



-^{E'f<^ + (^'A - i2n' + n'))'! - 2^nA'pf 
—{E'f^ - 2 [n" + nn') / + 2 (a - nA'f 
-^(h')2$ - 2 {n" + nn') / + 2 {2n' + n^) nf 

2A 



2A 



(h - 2 [n" - nn' - w) / 



(3.43) 



times 5jk- Here we have used (j3.14p twice and the definition $ := — (/' + nf). 

On the other hand, if we write ip = — nf on the left hand side of (|3.42p then we can replace V' by ^ 
plus the correction term 

- 2 {2n{nf)' + {nf)") = -2 {2n{nf)' + n"f + 2n'f' + nf) 

= -2 {n{nf)' + n"f + 2n'f' + n{f' + nf)') 

= 2n^' -2{nn'f + n"f + {2n' + n^)f') 

= 2 [n^' + {2n' + n^)^) - 2 [nn' + n"f - {2n' + n'^)n) / 

= 2 {n^' + {2n' + w^)^) - 2 {n"f - n'n - n^)n) / (3.44) 

times 5jk- We have again used the definition "3) := — (/' + nf) frequently. Combining (|3.43p and ()3.44p . we 
see that we can rewrite (I3.42P in the form 



[An^' + 2^" + 2{n^' + {2n' + n^)^) - 2{n"f - n'n - n^)f] 6. 



'jk 



.J1^0Z'-—v'{E)Z + C^ 
2A^ 2 ^ ^ ^ 



2{n" - nn' - n^)f 



Ojk 



(3.45) 
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or 



[27^(2^- + + 2^" + 2{2n' + n^)^] 6jk 




—v'{E)Z + iE')H 6,,. (3.46) 



Now the scalar contribution to (j3.30p is evidently given by 



= [.^-qi]^^^-[A{^-^) + 2{^" + {2n' + n^)<^ + ni2^ + <^y]5ab 

-^\2C'{<^ + Z')-v'(C)a''z]5ab. 



(3.47) 



This equation, of course, breaks into the two independent relations 



$ - * = 



(3.48) 



and 



[2^" + 2{2n' + n^)'^ + n{2^ + '^y]5ab = 2C'[<^ + Z']-v'{C)a^Z 6ab. 



(3.49) 



They are obviously equivalent to (|3.41|) and (j3.46p under the identification of the obvious variables and 
= H. Notice again how the conservation equation (|3.38p was used in order to establish this result. 

4 Comparison with SCPT Coupled to Dust 

In the previous section we showed that our manifestly gauge-invariant formalism with dust reproduces the 
usual formalism with linearly invariant quantities without dust, plus small corrections. In this section we 
construct, on the one hand, linearly gauge-invariant quantities with dust a la SCPT which are certain linear 
combinations of the dust, scalar and geometry perturbations. On the other hand, we expand our manifestly 
gauge-invariant variables, which are certain non- linear aggregates made out of the non- invariant dust, scalar 
and geometry perturbations, to linear order. We will show that the sets of these two quantities coincide 
precisely. 

The original action that we started with in the companion paper [T] is the action of the physical sys- 
tem "gravity -|- scalar field -|- dust". Then we rewrote the Hamiltonian and diffeomorphism constraints 
in an equivalent form such that the Hamiltonian constraint could be deparametrised and constructed the 
manifestly gauge- invariant observables Qij,P^^ and H,n, corresponding to the gauge- variant three metric 
Qab, the scalar field ^ and the conjugated momenta thereof. As a second step, the equations of motion for 
Qij,P^^ and H,n were perturbed around a flat FRW background, resulting in second order time derivative 
equations of motion for the perturbations of the metric, denoted by 5Qij, and the scalar field (5H. 
We discussed in the last section that, although starting with a system including the dust (and thus four 
more configuration degrees of freedom), the corresponding equation of motion for the physical degrees of 
freedom reproduce the results obtained by SCPT. As the latter does not employ dust, it deals only with 3 
physical degrees of freedom. 

In this section we want to consider the system "gravity + scalar field + dust" as our starting point again, 
but then follow the lines along the SCPT analysis and show that we can also use the perturbations of the 
dust fields, denoted by 6T,6S^, in order to construct an (up to linear order) gauge-invariant version of 
the perturbed four metric Sg^i,. In this case, we have also 7 physical degrees of freedom in the system (6 
gravitational and 1 matter). As a first step, we will use the dust fields 5T and 6S^ in order to derive ten 
gauge-invariant components of Sg^u and one for the scalar field 6^. Throughout the section we will denote an 
(up to linear order) gauge-invariant extension of a perturbed quaritity 5f by 6f. Secondly, we will demon- 
strate that 6 physical degrees of freedom contained in Sg^i, and 6^ agree with our observables constructed 
in the previous sections, if these are expanded up to linear order in the non-invariant perturbations of the 
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configuration variables. 



As was done in the last section, we will denote an infinitesimal gauge transformatiord by i-^ + 
where u = {u^, u°'). Here u"" = u°j_ + S"'^u,b is a proper three vector and the scalar function u is the solution 
of the equation u^ah^^"" = W^a- recall the transformation behavior of the single components of the 

perturbed four metric Sg^^, under such an infinitesimal gauge transformation. Throughout the section we 
will indicate background quantities with a bar such that 5/ := / — / for any quantity /. For the moment 
we will not decompose 5gfj,p into scalar, vector and tensor modes, because the whole analysis done in this 
section can be applied to the undecomposed tensor quantities. Denoting the gauge transformed quantities 
by ^5^1., we have 

(^500 = '^500 + {C-ug)m = Sgoo - 2a{au°y 

SgOa = SgOa + {Cug)oa = Sgoa + a^ (^-U°a + {u±+U'''ydabj 

Sgab = Sgab + {Cug)ab = Sgab + [l^U^bab + (^^1 + U'^) + {u\ + W^) ,^ba^ , (4.1) 

where the /' denotes a derivative with respect to conformal time, denoted by xq and defined through 



4.1 Review of the SCPT Strategy for "Gravity + Scalar Field" 

Before proceeding to the case of physical system plus dust, for the benefit of the reader let us briefly recall 
how the (up to linear order) gauge-invariant quantities are constructed in SCPT, in the case where no 
dust is present. A detailed revision of this construction can be found in appendix A of this paper. As 
a first step, one decompose the metric into scalar components denoted by (j),ip,B,E, vector components 
Fa, Sa and tensor components denoted by hab, where the latter is not of much concern, since it is already 
gauge-invariant. Their corresponding perturbed line element is then given by 

Sds^ = bg^^dx^'dx" = c? \^4>dT]^ + 2{B^a + Sa)dx''dr] + (2<^bab + 2E^ah) + Fa,b + F^,a + hab)dx''dx^ 

(4.2) 

from which we can read off the scalar, vector and tensor components of Considering how 5gf^i, 

transforms, we can derive the transformation behavior of 0, -0, E, B, Fa, Sa and (for completeness) hab- It is 
given by 

4)-^ = ^-uu^-[u^)\ ^'^ = ^ + nu^ 

E"" = E + u, B"" = B + u' 

5^ = Sa + iu'Jdab, F^ = Fa + u'i6ab 

Kb = Kb (4.3) 

where we introduced TC := a' /a. The strategy to construct gauge-invariant quantities in SCPT is (i) take 
the combination B — E' , which undergoes a shift in —u^ under a gauge transformation. This combination 
can then be used to construct gauge-invariant quantities for (p and ^. Recalling that we denote the (up to 
linear order) gauge-invariant version of any quantity / by /, we obtain 

4>:=^ = (j)-niB-E')-{B-E'y and := = + n{B - E'). (4.4) 



^Note, that in [S] the letter ^ is used for denoting their infinitesimal gauge transformation. Since ^ is already taken by the 
gauge-variant version of our scalar field, we decided to use the letter u'^ here, which should not be confused with the dust 
four- velocity given by capital V^. 
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(ii) One constructs a gauge-invariant vector by using the combination 

Va := Va = Sa- F',. (4.5) 

Consequently, we see that for the construction of gauge-invariant quantities in SCPT it is convenient to 
decompose Sg^^ into scalar, vector and tensor parts, such that four unphysical components of 6gfj,y^ namely 
E,B,Sa can be used to make the remaining six gauge- invarianllf]. For the benefit of the reader we also 
mention the capital letters $,^',14, used in the last section and in [9]. However, in this section we will 
use the "notation for (up to linear order) gauge-invariant quantities instead, because, when considering the 
gauge invariant quantity of 5g^y^ we prefer not to use 5G^y to avoid possible confusion with the ordinary 
Einstein tensor. 

4.2 Analysis for "Gravity + Scalar field + Dust" in Analogy with SCPT 

Let us now consider the system "gravity + scalar field + dust", for which we have 4 additional degrees of 
freedom, compared to the system "gravity + scalar field". In this case it is very natural to first keep all 
components of the perturbed metric 6g^i, and construct their gauge-invariant extension not by means of 
certain scalar /vector components of the four metric, but taking the four dust fields {5T,5S^) to perform 
this task. Since we want to derive Sg^y, there is also no need to consider the decomposed expression of 
that tensor. Therefore, after having recalled SCPT for the case of no dust, we will go back to equation 
(j4.ip . where the gauge-transformed components of 5g^u are displayed. In the following, we will show that 
the behavior of {5T,5S^) under infinitesimal gauge transformations suggests in a natural way that these 
fields can be used in order to make Sg^u gauge- invariant. For this purpose, we first need an Ansatz for the 

background dust fields (T, S^) compatible with the FRW symmetries. This can be derived from an Ansatz 
of the dust energy momentum tensor T^^,, specialised to the case of an FRW background. Recall that the 
general energy-momentum tensor was given by 

Tj^^^pU^U, with U^^-T^+(^^)S^;^, (4.6) 

where « indicates weak equivalence. That means that this equations is only true up to second-class con- 
straints. Ufj, denotes the four- velocity of the dust, while Pj, P are the momenta conjugated to T, , respec- 
tively. Note that here = (i, x"), since we have not yet transformed to conformal time . 
By considering the most general Ansatz for T^^ that respects the FRW symmetries, we get conditions on 
the dust four- velocity which carry over to the dust fields T, . The most general Ansatz is given by 

T'°=/o°Wa = and tS, = p°C7aC7b oc (4.7) 

In order to satisfy these requirements, we make an Ansatz for the dust fields {T,S^). For T := T{t) we 
take a, so far, arbitrary function of t so that T^a = 0. For = S^{x,y,z) we assume that it does not 
depend on time, so that 5^^ = 0. Furthermore, we require that 5 „ respects the homogeneity and isotropy 
requirements, that means S ^f, = 0. Therefore S = M^x" can only be linear in the spatial coordinates x", 
where is a matrix with constant entries. From the dust equation of motion shown in [1] in section 2, we 
know that CjjS ^ = 0. Inserting the Ansatz for T, S yields U"" = l/aP'Ua = 0. On the other hand, we have 
Ua '■= —T^a — ^S^a — ~^'^^a- One of our basic assumptions is that the matrix s'^^ with A;, a = 1,2,3 is 
invertible and we denoted its inverse by S^. Requiring now that Ua = —^S^a ^ and multiplying by the 

— a ~P 

inverse Sj on both sides yields the condition — =- = for all j = 1, 2, 3. Hence, the conjugated momentum 
Pk of vanishes in the case of FRW if T carries all the time dependence and 5^ all the spatial dependence 
^Recall that the vectors Sa,Fa are assumed to be divergenceless and therefore carry only two degrees of freedom. 
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of U^. Furthermore, from the fact that T defines proper time along the dust lines, i.e. C^T = 1 as a 
consequence of the dust equation of motion, we get (T^j)^ = 1, that is T^t = 1- Hence, we obtain T = t. 
Thus, the only non-vanishing component of is the it-component 

TS = p°, tI = tI = 0. (4.8) 

Thus, the symmetry requirements for tJ^^ lead to the usual four-velocity of pressure-free dust given by 
I7^ = (-1,0,0,0). 

Now, we want to analyse how the perturbed dust fields {6T,6S''), defined through 5T := T — T and 
5S^ := — S , behave under a gauge transformation x^^ — > x'^ -|- u^. The quantities {6T,6S^) will, 
for instance, be included in the perturbed dust four-velocity (5C/^ and thus in the perturbed energy mo- 
mentum tensoilll T^. For this purpose we will work , similarly to SCPT, in conformal time dx'^ = a^^dt. 
Then we have T ^ = aT^t = cl- Hence, the perturbations are given by 5T = T(t{x^)) — T{t{x^)) and 
6S^ = — = — M^x"". A gauge transformation u acts on 5T as 

= 5T + CyT = 6T + u^Tfl = 6T + u^t' = ST + au°, (4.9) 

where ' denotes derivative with respect to conformal time x^ as introduced already above. We can see that 
a~^5T is simply shifted by u^. In the last step we used that T = T q = a. in P] the authors use the 
combination B — E' which undergoes a shift of —vP under gauge transformations. Hence, in order to use the 
perturbations of the dust field, we can simply take a~^6T which is able to compensate any shift of in 6g'^^, 
in equation (j4.ip . The other shift which occurs for the components dg'^^ and 6g^^ are certain derivatives of 
the three vector u". This shift can be compensated by the transformation behaviour of 55^ which is shown 
below 

5(S^J" = 55,^, + {£uS'),a = SS':, + + v'X, = [5S^ + u%\a. (4.10) 

jf. , 

Here we used that S has no conformal time dependence. Consequently, the 5S themselves transform as 

5{S^Y = 5S^ + u^s'^t, = 6S'' + (n'^ + S^%,c)S% (4.11) 

which is the usual transformation rule for scalars under such a gauge transformation. Furthermore, we used 
the decomposition of u'' into a longitudinal and transversal part mentioned before. From the equation above 
we can read off that the combination 6S''s\.dab will be shifted by the vector Ua after having applied a gauge 
transformation. Therefore, this combination can be used to compensate for all shifts in Ua and derivatives 
thereof which occur in the gauge-transformed components of 6goa and 5 gab- Hence, the perturbed dust 
fields {6T,6S'') suggest very naturally a way to construct the (up to linear order) gauge-invariant metric 

perturbations Sgij,,y, since {a~^6T, s\6S^6ab) exactly provide the shifts in and u"", respectively, which are 
needed. The components of 5gf^u are, in particular, given by 

^goo = Sgoo + 2a{5T)fi = 5goo + 2a((5T)' 

Sgoa = Sgoa + a' [a~\dT),a - (sl5S''),oSab] = dgoa + [a~\dT),a - (sUs^Ydab] 

6gab = 6gab-a^[2^ST6ab + (Sl5S''),ahc + (Sl6S^),bSac]. (4.12) 

Note that all ten components of Sg^u are gauge-invariant now. These gauge-invariant combinations are 
different from the usual SCPT analysis where, as mentioned above, 4 out of the ten metric components 
are used to make the remaining six components gauge- invariant. Upon including the dust, the analysis 

^Note that, due to the vanishing of the background term Pk/P, the expression SS'^^j^ occurs in 5U^ only formally, since it is 
multiplied by the vanishing background quantity Pk/P- 
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above shows that the dust is a natural candidate for the unphysical degrees of freedom, because their 
changes under gauge transformation can in a simple way compensate for the changes of the components of 
6gfj_,y. A consequence is that when decomposing 5g^^ into scalar, vector and tensor modes, we also obtain 
gauge-invariant extensions of the quantities E, B, Sa, Fa- For completeness and to facilitate comparison with 
the standard SCPT framework, we also list the gauge-invariant extensions of these decomposed quantities 
E, B, Sa, Fa for the case where the dust fields are used to construct those gauge-invariant quantities. 
For this purpose, we need to decompose the vector 6w^ := Sl.6S^ into its longitudinal and transversal part, 
thus 5w'^ = 6w'^^ +(5?i;^ = 6'^'^6w\\^c + ^W^, where 6w\\ is the solution of the equation A(5w|| = 6w'^^. Explicitly, 
we have 



E 

Sa 



E 

Sa 



B 



^ = 11^- na-^5T 
--B + a-^5T-{5wu)' 



(t) + na-^5T + {a-^5T)' 

- {SvJ'JSab 

Finally, in order to derive the gauge-invariant extension of the scalar field perturbation 5^, we consider first 
the energy momentum tensor of the unperturbed scalar field ^, given by 

1 



Fa = Fa- 5wUab. 



(4.13) 



^rpmsitter 



(4.14) 



Here v{^) denotes the potential of the scalar field. Specialising to an FRW background as discussed in 
appendix A, we obtain for the components of T^^ 


1 



[matter 
00 

imatter 
Oa 



=matter 
^ ah 



;[(O' + aM0] 



^2— matter 



^2-matter^ 



(4.15) 



where we used p^^^^'^^ = 

7=matter 

1 components are 



Ui-^ + vii)) andp 



matter 



i(^2 -v{i)). The hnear perturbation jr™'^**'^'' of the 



rrpmatter 



jjimatter 



'-Oa 



rrrimatter 



Y \>^pSgab + ISabi^itfSgoo + 2tS^' - 



(4.16) 



Here we introduced v 



dv/d^. The gauge-invariant extension of 6T^^^^^^ , denoted by 6T^^^^^^ , can then 
be constructed by replacing the quantities ^^oo^'^^Oa in equation (j4.12p and 5^ by their gauge- invariant 
extension. We have not worked out the latter yet, so we need to consider how transforms under a gauge 
transformation. This is simply given by 



Consequently, the (up to linear order) gauge-invariant extension of 5^ is simply 



(4.17) 



6^ = 6^ 



-1 



Thus, the components of 6T. 



matter 



are 



rrpmatter 


1 

2A 


r/T^matter 
(>Toa 


1 - 

A . 


r/Timattor 
^Tab 


1 - 
A . 



^ S^,a + >^pSgOa 

Xp6gab + babi^itfSgoo + 2^5^ - a\-Jm] 



(4.18) 



(4.19) 
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This finishes our task to show that, by means of the four (perturbed) dust fields 6T,6S^, (up to linear) 
gauge-invariant extensions of the four metric Sq/j^u, as well as of the scalar field energy momentum tensor 
^jvmatter ^^^^^ Constructed. 

4.3 Comparison between SCPT + Dust and our Relational Framework 

As a last step, we will demonstrate that the (up to linear order) gauge-invariant quantities 5g^i, and ST^^^^^'^ 
derived in the last section agree with our (manifestly) gauge-invariant quantities 6Qij (+ corresponding lapse 
6N and shift 6Nj) and T™'^**™(5(5, 6'E), when these are expanded up to linear order. By this we mean that 
we will expand the formula in pLj for the gauge- invariant quantities F := Of, associated with any quantity /, 
around a given phase point (background), but now in terms of the non-invariant variables. The perturbations 
are then defined through 50 j := Of — Oj = Of — f, where we assume that the perturbation is around the 

configuration T = T, = . We consider only terms in which the perturbations occur at most linearly. 
Recall that in order to construct the (manifestly) gauge-invariant extensions of the perturbed three metric 
5qab, perturbed lapse function 6n and perturbed shift vector 5na, we first rewrote the Hamiltonian and 
diffeomorphism constraint c*°* = c^'^° + c^^^^^"^ + c° = c + and c*°* = cT° + c^^^^*^^ + = Cq + c^, 
respectively, in the following equivalent form: 



-tot ^ p _^ ^ ^-^j^ ^ _ Vc^- g""CaCb 



-tot 



C 



'3 



Pj + hj with hj := S^[-hTa + Ca\. (4.20) 



Note that the here has nothing to do with gauge invariance. We just took the notation used before in 
[1]. Since these constraints are mutually Poisson commuting, we could then perform our construction of 

the gauge-invariant extensions in two steps. First we construct diffeomorphism -invariant quantities O^^^ 
and afterwards use those in order to construct quantities that are also gauge-invariant with respect to the 

(2) 

Hamiltonian constraint denoted by F := O . The corresponding perturbed quantity is then given by 6F := 

(5(0^^(\)). Here, we are only interested in the terms of 5F which are linear in the perturbations, because we 

want to compare them with (up to linear order) gauge- invariant components of Sg^j^y in equation ()4.12p . For 
this reason, the general expression for 5F simplifies a lot and can be easily calculated. For the perturbations 
of the difFeomorphism-invariant observables 50fj we find up to linear order in the perturbations 



50f = sf+ I dMs''-s'){y){cl°Hy),f(.x)} 



6f- / dW(y){cr(y)>/(^)}- (4-21) 



Here we used 6f = f — f and, in notation introduced in [Ij, s'^ = o'^, as well as the definition of 6S^ = S'^ — S'^ 
in the last line. Note that up to linear order only Poisson brackets evaluated on background quantities occur, 
which we indicated by the bar above the Poisson bracket. This is due to the reason that 6S^ is already 
linear in the perturbations. Using 50^f\ we are now able to construct 6F, which up to linear terms in the 
perturbation given by 



6F 



50f+ I d^y{T-T){y){c'°'{y),60f\x)} 

6f- [ dW(y){ct°t(y),/(x)}- / d'y{T-T){y){c^-^{y),f{x)} 



6f - / d'y6T{y){c'-^{y),f{x)} - / d'y6S''{y){'^,-'{y), f{x)}. (4.22 
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Here we denoted [(5-F]iin, expanded up to linear order, by 6f as we did for 6g^,y, in order to emphasize that, 
in contrast to 6F, the former are not gauge-invariant up to ah orders. Furthermore we used that T = r in 
our case and 6T = T — T. We want to apply the above formula to / = qab-,^-, respectively, both of which 
are Poisson commuting with P, P^, the momenta conjugate to T and , respectively. Recalling the form of 
the constraints in equation ()4.20p . we see that equation (j4.22p simplifies even more to 

5} = 6f- I d^y6T{y){h{y)Jix)} - J d^y5S\y)Sl{ca{y)J[x)} (4.23) 

for those /'s which do not depend on the dust variables (T, S'^). Furthermore, we used that 

{{hTa){y),f[x)} = Ta[y){h{y)J[x)} = T,a{y){h{y), f{x)} = 0, (4.24) 
because we have T „ = 0. The computation of the remaining Poisson bracket above yields 

d^y6Tiy){h{y)Jix)} = f d'y6Tiy)(niy){ciy), f{x)} +na{y){cbiy), f{x)} 



dS5T{y)(niy){ciy),f{x)}) + (£^/)(x 



d-'y6S'{y)Sl{ca{y)J{x)} = (£^.55./)(x), (4.25) 

where we used the definition of the lapse function n := c/h and the shift vector := —Ca/h. Hence, the 
final version of the formula for 6f has the form 

6f = ^f- J d%ST{y) (n(y){c(y),/(x)}) - JcJ) - (^^^.Z). (4.26) 

In order to compare the results obtained from equation (j4.26p with the one for 6g^u derived in the last 
section shown in equation (j4.12p . we need to express the perturbations of the four metric in terms of the 
corresponding perturbations of the three metric Qab, lapse function n and shift vector n^. Recall that 
500 = o^(— + q°''^nani,), gQa = aua and Qab = qab- Consequently, the relation for the (up to linear order) 
gauge-invariant perturbations is given by 

S^jm = 2n(5n — 6qabn°'n^ + '2^^na5nb) = —2c?Sn 

5goa = a5na 

Sgab = Sqab, (4.27) 

where we used that n = 1 and n"" = 0. Starting with the spatial - spatial component of the metric Sqab and 
applying the formula in equation (j4.26p to it, we obtain 



lab 



Sqab - j d^y5T{y)[n{y){c{y),qab{x)]^ - (^n9)„f, - (^55*9)^^ 

= Sqab - j d^y6T{y) (^{c{y),qab{x)}^ - JC^^ (4.28) 

Here we used in the last line that for FRW the lapse function and the shift vector are given by 

n = cjh = cj \J (? — q°-hcaCb = c/c = 1 and Ha = —Ca/h = —Caj \/ (? — q°^CaCb = 0. (4.29) 
The remaining integral in equation (I4.28[) yields 



/ 



d^ydT{y)(^{ciy),qabix)}') = j d^y5T{y)[-^={qacqbd -\qabqcd)p'"'){x)5{x,y) (4.30) 



5T{x)^{-]^a^){-2a)5ab 

6T{x)2aa 

2a'dT{x). 
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In the second line we inserted the exphcit form of for FRW and the corresponding momentum p"'^ 
conjugate to it, given by = a'^6ab and p"^ = —2a6"'^, respectively. In the last step we replaced the 
derivative with respect to T by one with respect to conformal time. Considering the explicit form of 
also for the Lie derivative occurring in equation (j4.28p . we end up with 

Sqab = Sqab-2a'6T6ab-a^(Sl5S''),aScb-a\Sl6S''),bSac 

= Sqab - [2^6T5ab + (Sl5S^),aScb + (Sl5S^),bSac] ■ (4.31) 

The last line agrees exactly with the expression in equation ()4.12p if we identify 6 gab with 5qab- Next, we 
want to discuss the case for the scalar field perturbation 6^,. Since any spatial derivative of the background 
scalar field ^ vanishes, the Lie derivative with respect to Sf^dS'' in equation (|4.26|) . evaluated on the FRW 
background, vanishes as well. Therefore the (up to linear order) gauge-invariant extension is simply 

6^ = 5i- j dVny)({c(y),/(x)}) (4.32) 



5i- j d'y5T{y)(^-^'K){x)5[x,y) 



1 



= 5i-5T{x)[^==i^\ 
^ Vdet(g) 

= 5i- a~^t^T{x), (4.33) 

which exactly coincides with the result in equation (|4.18p . Note that we reexpressed the momentum vf 
conjugated to ^ in terms of ^ in the third line, which in the case of FRW has the simply form vf = ^Vdetf. 



Finally we have to discuss the lapse function 5n and the shift vector 5n°. We want to stress that in 
the rest of these two papers we always have treated them as pure gauge with respect to the primary con- 
straints and have reduced the problem to the remaining degrees of freedom. They never appeared anymore 
as independent variables and in the effective, physical Hamiltonian description lapse and shift were derived 
quantities as well, they were not independent. 

4.4 Subtleties with Lapse and Shift Functions 

In order to compare our framework with SCPT we must deal with the following subtlety: 

As mentioned earlier, in the Hamiltonian framework we have ten gravitational configuration and momentum 
degrees of freedom to begin with. These are n, n", qab, denoting lapse function, shift vector and three metric, 
respectively, as well as their conjugate momenta p, pa, respectively. This phase space is subject to four 
primary constraints z = p = 0, = = 0, as well as four secondary constraints c*°* = c'^^^^ +0^*^° _j_gmatter 
and c*°* = c^^^ -\- cf^° -\- q^^^^^v ^ These eight constraints are first-class and play a dual role: on the one 
hand they constrain the phase space to the constraint surface, on the other hand they generate gauge trans- 
formations via canonical transformations. The constraint equations can be used to eliminate eight of the 
momenta while the gauge transformations eliminate eight of the configuration degrees of freedom. Hence 
the physical or reduced phase space has eight canonical pairs less than the unreduced one. The primary 
constraints and the gauge transformations generated by them, which consist in arbitrary changes of lapse 
and shift, completely decouple from the rest of the equations, namely c*°*, c*°* do not depend on n, n"',p,pa, 
in particular, the gauge transformations generated by c*°*, c*°* do not affect lapse and shift functions. It is 
therefore mathematically convenient to simply forget about p,Pa and to treat n,n'^ as Lagrange multipliers 
which is customary. 
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Let us compare this with the situation in the Lagrangian framework. Here we consider only the ten grav- 
itational configuration coordinates, momenta are never introduced. The system is subject to Einstein's 
equations. Four of them, the temporal - temporal and the temporal - spatial ones are equivalent to the sec- 
ondary constraints that we also find in the Hamiltonian formalism. Namely, in the Hamiltonian formalism 
the equations of motion are generated by the primary Hamiltonian 

i^primary = / X {l^ Z + 1^'' Za + n C*°* + n'^ 0*°*) , (4.34) 

Jx 

where v, v'^ are Lagrange multiplier fields. The equation of motion for qab following from ()4.34p enables us 
to express in terms of qah, Qab, n, n", specifically 

p'^* = v/dit(^ (g - g'V'') i^cd with Kab = ^[qab-2D(,nb)). (4.35) 

Here we introduced the extrinsic curvature Kab, and D is the torsion free covariant differential compatible 
with Qab- When inserting (|4.35p into c^^°, ci'^°, one obtains precisely the temporal - temporal and temporal 
- spatial components of the Einstein tensor, respectively. Thus we can match four of the equations of motion 
in the Hamiltonian and the Lagrangian equations of motion. The remaining spatial - spatial components of 
the Einstein equations are obtained in the Hamiltonian formalism by inserting (j4.35p into the equations of 
motion for (jab which also uses the equation of motion for p"'' as obtained from ()4.34p |19j , so that we get an 
exact match at the level of the equations of motion (one has to use that the spacetime metric is expressed 
in terms of the configuration coordinates as gu = — -|- qabn""^^, 9ta = Qabi^^, 9ab = Qab)- Now in the 
Lagrangian formalism the secondary constraints do not generate gauge transformations because one does 
not have a phase space formulation. Rather one notices that the Lagrangian and thus the Einstein equations 
are invariant under spacetime diffeomorphisms. Thus one imposes spacetime diffeomorphism invariance as 
a gauge symmetry. Now we have seen that the temporal - temporal and temporal - spatial components 
of the Einstein equations do not involve time derivatives of lapse and shift. Thus it is natural to eliminate 
them via the secondary constraint in terms of qab, Qab and to eliminate four components of qab by using 
spacetime diffeomorphism gauge transformations. 

We see that in both formalisms we get the same number of physical configuration degrees of freedom, 
however, the avenue to get there is somewhat different. In particular, there is no counterpart of the gauge 
symmetry associated with the primary constraints at the Lagrangian level. Furthermore, while the secondary 
constraints generate spacetime diffeomorphisms on qab [H], they preserve lapse and shift fields, whereas in the 
Lagrangian formalism lapse and shift do change under spacetime diffeomorphisms. Thus, there is a mismatch 
in the number and action of the gauge transformations of the Lagrangian and the Hamiltonian formalism 
at the level of the unreduced configuration space. It is clear that the physical degrees of freedom that 
both formalisms produce are the same, since both formalisms tell us that one should extract the spacetime 
diffeomorphism invariant information contained in qab (as far as the gravitational degrees of freedom are 
concerned). However, before using the secondary constraints, the diffeomorphism invariant quantities that 
one constructs in the Hamiltonian theory do not involve lapse and shift because they are invariant under the 
secondary constraints, while in the Lagrangian formalism one has to use them. This has the consequence 
that in the Lagrangian formalism one has in principle six rather than two physical configuration degrees 
of freedom whereas in the Hamiltonian formalism one has only two. In the Lagrangian formalism, four of 
them are eliminated via the secondary constraints and thus, at the end of the day, one gets an exact match 
at the level of the physical degrees of freedom. But before doing that it is somewhat difficult to compare 
the gauge invariant degrees of freedom in the two formalisms. The reason for why that happens lies deeper 
and has to do with the question whether Lagrangian (Noether) symmetries have canonical generators in the 
Hamiltonian formalism. This is analysed in all detail in [20j. 

Fortunately, there is an elegant solution. The idea |20] is to add, in the Hamiltonian formalism, to the 
secondary constraints a linear combination of the primary constraints with carefully chosen coefficients such 
that the new constraints generate spacetime diffeomorphisms on the full phase space, including lapse and 
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shift. This leads to an equivalent set of constraints which leads to the same final set of physical degrees 
of freedom but which has, after reducing with respect to the new secondary constraints only, the same 
gauge-invariant degrees of freedom as in the Lagrangian formalism. Then afterwards one must impose the 
secondary constraints in the Lagrangian formalism and reduce with respect to the primary constraints in 
the Hamiltonian formalism in order to arrive at the true degrees of freedom. We will now sketch how this 
is done as we will need this in the subsequent subsections. 

Fortunately, we do not need all the machinery developed in [20], but can rather develop the required 
formulae by basic methods. The starting point is to recall [19] that for arbitrary test functions b, 6" we have 

{^"\b) + ^°\b),q^^{x)}=Cuq^,u{x) where u^" = bn^" + X^^b" (4.36) 

with c*°*(/) = d?x b{x) c*°*(x) and c*°*(6) = J-^ (Px b°'{x) c|f*(x). Here Cu denotes the spacetime Lie 
derivative with respect to the vector field u and q^u = g^y + n^rii, is the three metric on the leaves Xt of 
the foliation, with unit normal n'^, parametrised by the one parameter family of embeddings Yt = Y{t, .) : 
X ^ Xt- In deriving (j4.36p the equations of motion have been used and thus we can say that the secondary 
constraints generate infinitesimal spacetime diffeomorphisms parametrised by u. Notice that any vector field 
can be split into components tangential and orthogonal to the foliation so that this is no loss of generality. 

We now would like to add to c*°*(6) and c*°*(6), respectively, a term linear in p and pa-, respectively, 
with smearing functions as coefficients that are linear in both 5, b"" such the extended constraints also 
generate infinitesimal spacetime diffeomorphisms on lapse and shift functions parametrised by the above 
vector field u. In order to determine those coefficients, we first need to determine the transformation 
behaviour of lapse and shift functions under infinitesimal spacetime diffeomorphisms. Starting from the 
identities gu = —n^ + qab^'^n^, Qta = lah^^^ 9ab = Qab^ we can solve for lapse and shift in terms of spacetime 
metric components 

= g'^'gtb, = -9tt + g^^gtaQtb, (4.37) 

where g""^ is the inverse of the spatial metric gab ■ Now using that the change of the metric under infinitesimal 
diffeomorphisms is given by the Lie derivative 

KOiiu = {^vg)/,u = uPgi,u,p + 2gpi^^uP^y (4.38) 

we can work out the infinitesimal transformation of lapse and shift. The result is, in the frame where X^ = t 
and X"" = x"" (frame adapted to the embedding), 

= u^rff, - n^uffe + nf^ + (n"n*),t - (n^g^^ + n"n'') u*^ 
InbuU = -n"-n^ {vfgab,c + 2u'^(^^gb)c + u^gab,t + 2nan*b^ 

+2n" (u^gab,t + U^^tdat + gttu\a + y^gatfi + gabU^,t + u'^^agbc'Tl''^ 

- {gtt,tu^ + 2gttu\t + gtt,a + ^Uau'^t) ■ (4-39) 

Here we have used the explicit form of ()4.38p . Formula ()4.39p is rather complicated. Worse than that, it 
contains time derivatives of n, n". We cannot use the equations of motion in order to replace those time 
derivatives by the momenta p, pa because the very reason for the appearance of the primary constraints 
p = Pa is that the Legendre transform is singular, as the Einstein Hilbert action does not contain time 
derivatives of lapse and shift. This is different for time derivatives of qab- Hence, in order to write (I4.39|] 
in terms of canonical coordinates, we must get rid of the time derivatives of lapse and shift functions. 
Miraculously, the equations ()4.39p simplify drastically when we split u into normal and tangential pieces, 
= bn^ + b^-y^a- Noticing that in the chosen frame we have from g^^n^n'^ = —1, Y^ag^u^'^ = that 
n^^^ = 1/n, n^^"" = —n"'/n, as well as nt = —n, Ua = 0, we find 

h n°- 

= -, u" = b'' -b—. (4.40) 

n n 
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Substituting (j4.40p into (|4.39|) . almost everything cancels due to the explicit lapse and shift dependence of 
the vector field u and yields after some calculus 



Sun'' 
2n6un 



(4.41) 



The time derivatives of lapse and shift have disappeared as desired. The required extension of c*°*'(6) + c*°*(6) 
that generates spacetime diffeomorphisms on all ten components of the spacetime metric and not only its 
spatial - spatial components is therefore given by 



X 



{C^nf + q'^'ibnb - nbb) + b'l +P [b,t + b'^Ua - n'^ba] 



(4.42) 



The coefficients appearing in (j4.42p are no accident, but directly derive from the Dirac algebra of secondary 
constraints as shown explicitly in [20j. It follows that the primary constraints have to be there in order to 
make the formalism manifestly spacetime covariant. We observe that when choosing b = n, b"" = n"" the 
modified expression c^°*(/) + c*°^(/)+ (l4.42|) turns precisely into the primary Hamiltonian (I4.34|) . because the 
equation of motion for lapse and shift can be used to eliminate v = b, v"" = W. Thus we may say that for those 
special smearing fields the generator of spacetime diffeomorphisms is nothing but the primary Hamiltonian. 
Thus, the extended generator of spacetime diffeomorphisms is nothing but the natural extension of the 
primary Hamiltonian to arbitrary test functions different from b = n, b"" = n"". 

Performing integrations by parts (using suitable fall-off properties of 6, 6"), we collect the coefficients 
of 6" and b in ()4.42j) (but not of b"", b), which give the lapse and shift (LS) contributions to the spatial 
diffeomorphism and Hamiltonian constrain, t respectively: 



Pb n 

j^„ab 



,a 



Pbn,a- 



(4.43) 



The first equation tells us that n, n'^ transform as scalar and vector, respectively, under spatial diffeomor- 
phims, while p,Pa are correspondingly scalar and covector densities, respectively. For this reason all the 
derivatives in the second equation in (|4.43p can be replaced by covariant ones. 

Since we have extended the phase space by lapse and shift, we now must include the contributions (I4.43P 
into the secondary constraints which thus read 



tot/ 



.dust 



matter! _, dust 
J ~r '-a 



=: c' + d^""^^ 

-■ ^' _|_ --dust 



(4.44) 



They have to be used in our formula for gauge-invariant completions with respect to the secondary spatial 
diffeomorphism and Hamiltonian constraints, respectively. Now, as before, we solve ()4.44p for the dust 
momenta and obtain 



cr = P^ + S';{-h'Ta + (^a) (4.45) 

which is completely analogous to the formalism in which the primary constraints have been reduced already. 
The b and 6" terms in (|4.42p cannot be taken care of in c^^,c^^. In principle, we have to consider 
:= b, v"" := 6" as independent test functions. Therefore, the bp + b^'Pa = i^'^z^ contributions in ()4.42p 
have to be attributed to the primary constraints = p^, that is, z = p, Za = Pa- Since the primary 
constraints have not yet been reduced, we must also perform the gauge-invariant completion with respect 
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to the primary constraints. Notice that since z, Za, c*°*, c*°* are weakly Poisson commuting, so are the 
z, Za, c*°*', 0*°*' because ()4.43p is a hnear combination of the z,Za- Hence the system remains first class. 
Since z, Za are presented in deparametrised form, we know that the constraints z, Zq, c*°*', cjf*' are mutually 
Poisson commuting. Hence, in addition to what we did in our companion paper, we must supplement the 
gauge-invariant extension F of some function / with respect to the secondary constraints by gauge-invariant 
extensions with respect to the primary constraints z, Za, for which of course we use as clocks the lapse and 
shift functions n^nP". Specifically, that last extension is given by the formula 



0f[v] = F + Y,-^ I d^xi.J d^xn K^-n^^)(xi)..K"-n^")(x„)..{z^,(xi),..,{z^„(x„),F}..} 

n=l -^-^ -^-^ 

(4.46) 

with arbitrary functions z^^, /x = 0, 1,2,3 and we have set := n, n'^^"' := n", o = 1,2,3. This operation 
leaves functions independent of n,n^ unaffected, of course, and, more generally, replaces any by ly^, in 
accordance with the fact that n, are pure gauge with respect to the primary constraints. 

Notice that the gauge-invariant projection O^^^ has to be performed after the operations 0^^\ O^^^ of 
our companion paper, because p,Pa Poisson commute with T,S^, but n, do not Poisson commute with 
cjf*', cjf*'. Since in the complete projection the constraints never act on the r — T, — , v — n, u"" — n°' 
powers, (j4.46p has to applied in this order. 

With this machinery at our disposal we can now continue our comparison with the SCPT framework. 
4.5 Comparison with the SCPT framework for Lapse and Shift 

As discussed in the last section, for the purpose of comparing with the SCPT framework, we have to use 
the extended Hamiltonian formalism. Thus we now must use the constraints shown in equation (j4.44p . as 
well as the additional projection in equation (I4.45p . All that we have said before remains valid because on 
the constraint surface defined by the primary constraints the constraints ()4.44p coincide with the old ones. 
Let us write the perturbed expansion to linear order of a complete gauge-invariant quantity in the extended 
phase space as 



SOf = 5f- / d'y i6T{y) {c'-^'{y), f} + 6S^ (y) {cf (y), /} + <5n(y) {z{y), f} + dn'^iy) {za{y),f} 

(4.47) 

where 

6T = T-T, 6S^ = S^-a^, 5n = n-v, and Sn" = n" - v"" . (4.48) 
Evaluating (j4.47p for the functions / = n(x), n"'{x) yields 

5n{x) := 50^^^) = \dn - ^6S^n^ + [6T]^^ - 6n\ (x) (4.49) 



q^b (5Tn^b-n[6T] 



6n^ 



(x) 



(4.50) 



where we used n = 1, n° = 0. We want to compare n and fia with Sgoo and 5goa, respectively, in equation 
(j4.12p . From equation (j4.27p we can read off 



6g 



00 



-2a^(5n and Sgoa = a^Sah^n^ 



(4.51) 
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Hence the expressions for J^qo aiid 5gQa in equation (j4.12p are equivalent to 



6h = 6n 



6T 



and = Jn" 



(4.52) 



where we used that 5^ 



0. Comparing the Lagrangian result (|4.52p with the Hamiltonian result (|4.49p . 
we do not seem to obtain a match. The reason is the appearance of the time derivatives of 5T and 5SK 
On the other hand, notice that with the identifications b = 5T and 6" = Sj5S^ the square bracket terms 
on the right hand side of ()4.52p precisely cancel the spacetime diffeomorphism transformation derived in 
equation ()4.4ip via Hamiltonian methods. The difference is that the fields b, b"" in equation (j4.4ip were 
smearing fields, while here 5T and 6Sj6S^ are phase space dependent functions. The only way in which the 
Hamiltonian formalism can produce time derivatives of the canonical fields is via the equations of motion. 
Now the equations of motion in the Hamiltonian formalism which lead to an exact match with the Euler - 
Lagrange equations are with respect to the primary Hamiltonian (I4.34|) . We find 



T'(x) — {-f^primary; ^(2^)} 



n"T„ - n 



P 



P2 + gabc'^c', 



S^{x) 



{^primary, T(x)} = [u^ S^^,] (x). 



Perturbing equation (I4.53P around the FRW background and using 

5f = 5n and 'S^5S^ = (5n". 



T 



0, as well as 5;^ 



(4.53) 

0, we find 
(4.54) 



Equation (j4.54p is evidently the missing link to bring (j4.49p and (j4.52p to an exact match. 

By construction, the full, non-perturbative expression for 0„a(^) and On{x) will be expressed in terms 
of the fundamental seven physical degrees of freedom corresponding to qab,C and their conjugate momenta. 
Hence in the Hamiltonian framework ()4.52p are of no further interest. In the Lagrangian formalism, the 
secondary constraints fulfill the task to express those equations in terms of the other degrees of freedom. 
Hence in both frameworks we end up with seven gauge-invariant degrees of freedom. 



Summarising, we showed that, using the four dust fields 6T,SS^, we are able to construct (up to linear 
order) gauge-invariant quantities along the lines of SCPT. We obtain 10 gauge-invariant components dgfj,i/ 
and one for the matter scalar field S(,. Out of these 11 degrees of freedom only 7 are physically relevant, 
because apart from the non-physical dust degrees of freedom the system is reduced by another four degrees 
of freedom due to the four primary constraints of General Relativity. This leads to a reduction from the 
15 configuration degrees of freedom of the system "gravity -|- scalar field -|- dust" down to 7 true degrees 
of freedom. These seven degrees of freedom reside in the perturbations of the three metric 6gab = ^Qab and 
the scalar field 6^. Finally, we expanded the (manifestly) gauge-invariant observables constructed by the 
method introduced in [Ij up to linear order in the perturbations and compared them with 5g^i, and 5^. It 
was shown that these two ways of constructing (up to linear order) gauge-invariant quantities lead to exactly 
the same results. Four of those variables are redundant in both formalisms. Thus, for the physical degrees of 
freedom we get an exact match. In contrast to section [3] where SCPT without dust was considered, we also 
have an agreement on the number of physical degrees of freedom here, 7 for SCPT and 7 in our framework. 
However, as was discussed in section [3l the equations of motion derived for the 7 degrees of freedom show 
that 4 of these degrees of freedom freeze out in the late universe, so that the physics is mainly governed by 
the 3 degrees of freedom used in SCPT. 



5 Summary and Conclusions 

In this paper we applied a new framework for general relativistic perturbation theory, developed in our 
companion paper [1], to the important case of cosmological perturbations. The central feature of our new 
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approach is the use of a dynamically coupled observer medium, given by pressureless dust. This allows for 
a complete deparametrisation of the physical system of interest, that is gravity coupled to whatever matter 
one wishes to include. As a result, the usual gauge freedom of general relativity is eliminated and true 
observables are obtained, together with a true time evolution generated by a physical Hamiltonian. By 
specializing our general framework to the case of an FRW background spacetime, we developed a manifestly 
gauge-invariant cosmological perturbation theory. As the quantities that are perturbed, namely the three- 
metric and all the (non-dust) matter fields, are already fully gauge-invariant by construction, the familiar 
problems with gauge freedom, that had been troublesome for standard cosmological perturbation theory 
(SCPT) for a long time, never arise here. In particular, it is straightforward, though no doubt involved, to 
derive perturbed equations of motions to arbitrarily high order. 

In this paper, we limited ourselves to investigating the linear order. We found that up to a small 
correction term, our formalism reproduces the known SCPT equations. The correction term is an imprint of 
the dust system. On physical grounds, this could have been anticipated right from the start, as the observer 
dust adds to the overall energy momentum balance. Our framework thus illustrates the fact that including 
a realistic, that is non-idealized, observer in the description of a physical system will always disturb the 
original system to some extent. As this correction turned out to be inversely proportional to the scale 
factor, however, it decays away very quickly in the early universe. Consequently, the equations of motion 
coming out of our framework are well within an acceptable range of the standard cosmological model results. 

There are many routes for further investigations that one can take from here. Since the issue of gauge 
invariance has been entirely settled to all orders in our approach to perturbation theory, the obvious next 
step is to calculate the equations of motion for higher order perturbations, in particular for the second order 
|18j . As mentioned in the introduction, the latter is a topic that currently attracts a considerably amount of 
attention, due to its connection with the issue of non-Gaussianity of cosmological perturbations. This should 
also help to settle the stability issue of linear perturbation theory |23J. Another interesting project is to 
investigate the spherically symmetric sector of the theory and to determine whether there are modifications 
of Newton's law at large distances. In fact, it might not be totally surprising if there were modifications 
because the dust is a perfect candidate for a WIMP, one of the candidates for dark matter, since it only 
interacts gravitationally. 
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A Review of Standard Cosmological Perturbation Theory 

For the benefit of the reader we collect here the most important formulae of the SCPT formalism, adapted to 
our notation. In particular, notice that we are using the relativist's signature (— , -|-, -|-, -|-), which is opposite 
to the one used by cosmologists. 

A.l Curvature and Energy Momentum Tensor 

Our convention for the curvature tensor is 
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Its explicit expression in terms of the Christoffel symbols (we use the notation (.)_^ := d{.)/dx'^) 

^iiv — 9 ^^pixvi ^pixv = 2 {9piJ.,v + 9pv,ti ~ 9iJ,v,p) 5 (■^■2) 

corresponding to the torsion free covariant differential V compatible with a general metric g^y with inverse 

R^..p " = 2 (-5[^r-]^ + rj[^ r-],) . (A.3) 

Prom this we infer the Ricci tensor, Ricci scalar and Einstein tensor, respectively: 

Rpiy '■= Rppv R '■= 9^'^ R^ivi Gfj_u = R^y — -Rg^y. (A. 4) 
The energy momentum tensor for bosonic matter is given by 

q-i 2 ^'S'mattcr ( \ K\ 

vm9T\ ^ ^ 

which, for a minimally coupled scalar field ^ with action 

Scatter = "7^ / V\d^\ {g^'' {V yQ + v{0) (A.6) 

(where v is its potential), becomes 

XT^y = [V^C] [V.C] - ^5^. (5''"(VpC)(V.C) + v{0) . (A.7) 
The Euler-Lagrange equations for the gravity and scalar matter system (including a cosmological term) 

5geo = - / d'^X v/|det(5)| {R - 2A) (A.8) 



are the Einstein equations 
and the matter field equation 



Gf,y + Ag^^y = '^T^y (A.9) 
g'^''V^VyC=^v'{C). (A.IO) 



A.2 FRW Background 

The spatial constant curvature k = FRW line element with scale factor a{t) reads 

ds'^ = -dt^ + aitfSahdx^dx^ = aiaPf-q^ydxi'dx'' (A.ll) 

where we have switched to conformal time dx^ := dij := dt/a{t) and r]ny denotes the Minkowski metric. In 
the second step we should actually use a different symbol d{x^) := a{t{x'^)), but we will slightly abuse the 
notation as is customary. 

All quantities that refer to a background metric will carry a bar. For the FRW background metric 
Vfiu — '^^ViJiv a short computation reveals 



i 00 — ~ ~- 

r° - 

Too = 

ro6 = ns^ 



rl = 0. (A.12) 
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Here and in what follows a prime denotes derivation with respect to conformal time x = rj, rather than 
cosmological time t, and we have introduced the Hubble function 7i = a' /a. It is useful to note that ()A.12p 
implies that 

Tl^ = An5^o. (A.13) 
Carefully using the definitions, we find again after a short calculation 



Roo = 


-sn' 


Roa = 





Rab = 


{n' + 2n^) 5ab 


R = 




Goo = 




Goa = 





Gab = 


- {2n' + ri^) 6 



The non-vanishing components of the energy momentum tensor are 



^Tab 



1 



1 



{Cf -a'v{C)]5ab=:a^\5abP, 



(A.15) 



where we have introduced background energy density p and background pressure p, respectively. The matter 
field equation becomes 



A. 3 Linear Perturbations 



1 



4(«V)' 



1 



«'(C). 



(A.16) 



We consider perturbations Sg^y := g^y — g^^ and 5C, = C, — Q. Any function F = F{g, Q of the metric g and 
the scalar field C is expanded to linear order in 5g, 5C,, that is, 5F will denote the linear order term in the 
Taylor expansion of F{g, Q) — F{'g, C). For instance. 



pflU 



1 



pu,fj, '~'ypu,p 



For the FRW background this yields 



5rz 



1 



-r]" SgxpT^^ + iff -{6gpp,u + 5gpu^^ - 5g^y^p) 



(A.17) 



(A.18) 
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Inserting (jA.lSp . we find explicitly 

00 



"J- ab 



1 

2^ 
1 

1 

2^ 
1 



[2n6goo6ab - {25go(a,b) - ^d'ab) 
{-2n6goa + 26g'Q^ - 6goo^a) 



S'^bO = ^ {-2n6gab + 5g'ab + '^^%[a,b]] 



^^bc = — [-T~(-^90aSbc + '^^9a{b,c) ~ ^9bc,a 



It is useful to note that (|A.19P implies 



^ [s^^Sgbc - Sgoo) 



Now we use (jA?2]) and fO]) to find 

6R 



(A.19) 



(A.20) 



(A.21) 



Specialising (|A.21[) to the temporal ~ temporal, temporal - spatial and spatial - spatial components, respec- 
tively, yields 
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5Rnh 
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1 



—,A5goo - n ( ^6g,,6^'' ] - SH [ —,6g, 



2a? 



1 



2a2 



yoo 



2a2 



2a2 



red i 



n 



1 

~2^ 

+2n 



S'"^5gcd - Sgi 



00 



+ 



^'^SOO'^afe 



-^^900Sab 
1 

2a2' 



1 



2a2 

^ 1 

2a2 



2a2 



+ 



(2(550(a,fe) - ^9ab) 



-2n ( -^<55afe ) + n6ab ( x^(<5'''<55cd - 5ffoo; 



(A.22) 



Here A = 5°'^dadh denotes the flat background Laplacian and 25Tcab = 2(5 (7^(0,6) ~ Sgab,c- This implies 



5R 



-r^r' 6g,^R.,+r''SR, 



1 



j2\xcd 



4 ( -3n'6goo + in' + 2n')5"'5g,d ) + 4 ( --5^00 + ^^''SRcd 



1 / \ f Ti 6"^ ^ 

+_ (-3W(5^'^(55oc,d + S'"'6''''6T,,b,d) + 2n U^Sgoo - ^(2<5<7o(e,d) - Sg',^) 



+3n { ^{6^%ed - Sgoo) 



(A.23) 
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Since 

we find 

6 Goo = 

SGoa = 
SGab = 



R 5 



(A.24) 



5Roa-^{rL' + rL') 6goa 



in' + n^)6'%cd + 3n^6goo^ 



(A.25) 



SRab -—in' + n^)dgab - —Sab 

2 



^4 ^ -3n'6goo + {n' + 2?^2) S'^dg.d) + ^ (-<5i?oo + <5'='''J^cd 



Next the perturbation of (IA.7P yields after a short computation 

AJT^, = 2C'6li5C),,) - \xp5g,,, + ^77^- {l^{^)Hg^^ + 2C'5C' - a'«'(C)^c) (A.26) 



or, in components, 



SToo 

SToa 
STab 



l.(2C'sC' + a\'i05C-v{C)5goo 

1 



A 



XpSgab + ^Sab (^-^{('fSgoo + 2C' 8C' - a\'(C)<5c) 



(A.27) 



A. 4 Parametrisation of perturbations and linear invariants 



The complicated equations of the previous section can be decomposed into scalar, vector and tensor con- 
tributions. As in [9j we introduce four scalar fields cj), ip, B, E, two transversal (with respect to the flat 3d 
Euclidean background metric) covector fields Sa, Fa and a symmetric tracefree, transversal tensor hab- Thus 
S"'^Fa,b = S°'^Sa,b = S^^hab,c = ^"'^hab = 0. These fields encode the ten independent components of dg^j.^ as 
follows: We write the perturbed line element 

ds^ = a^rjf.^dx^dx'' + 6goo{dx^f + 25goadx°dx'' + dgabdx^dx^ (A.28) 

in the form 



ds^ 



(-1 + 2(t)) {dxy + 2{Sa + B^a) dx^dx"" + ((1 + 2i;)6ab + 2E^ab + 2F(„,b) + hab) dx'^dx' (A.29) 



from which one reads off 



^500 = +2a 4> 

Sgoa = {Sa + B^a) 

SgOa = [2{^p5ab + E^ab + F(a,b)) + hab] ■ 



(A.30) 



Notice that we use different signs from the literature, due to our signature conventions: All perturbations 
enter with a positive coefficient proportional to a^. One speaks of scalar, vector and tensor perturbations, 
respectively, when only the respective field perturbations are non-vanishing. 

The fields (j),ij),E,B,Sa,Fa,hab are not invariant under (infinitesimal) gauge transformations, that is, 
spacetime diffeomorphisms. In general, a metric changes under an infinitesimal diffeomorphism generated 
by a spacetime vector field n'^ by the corresponding Lie derivative 



t=0 



(A.31) 
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where 1 1— > (/jj* is the one parameter family of diffeomorphisms generated by the integral curves of the vector 
field u [19j. If we take the order of u to be the same as the order of the general perturbations, then we see that 
we can gauge away four of the ten perturbation fields. We define := rj^^u'^ , that is uq = — u", Ua = u'^ 
and u := A^-'^n"^, = Ua — daU, where A~"^ is the Green function associated to A. We find explicitly 

SuQm = 2a^ ( - {au)' 
\a 

6ugOa = «^ ( [""0 + u] + Ua' 

SuQab = 2a^ {-HuQdab + u^ah + ufa,b)) ■ (A.32) 
Comparing with ()A.30p . we read off 

5u<p = - (ano)' 
a 

5uB = uq + u 
5uip = -Huq 
6uE = u 

SuFa = 

6Kb = 0. (A.33) 

From ()A.33P we immediately see that a complete and linearly independent set of linear invariants under 
spacetime diffeomorphisms is given by 

^ = ^-l[a{B-E')]' , ^ = ^ + n{B-E'), Va = Sa-F',, hab, (A.34) 
a 

which are six gauge invariant degrees of freedom. This follows from the fact that 

6uiB-E') = uo, 6uE = u, 5uFa = ui. (A.35) 

A scalar field transforms under diffeomorphisms as 

5uC = u^Cm, (A.36) 

which to linear order equals vPC^ = —C'uo- Thus the linearly invariant scalar field perturbation is simply 
given by 

Z :=6C + t{B -E'). (A.37) 

Just like the metric g^^, any symmetric tensor such as the Einstein tensor or the energy momentum 
tensor T^jy transforms as in ()A.3ip under infinitesimal spacetime diffeomorphisms. Therefore, the Einstein 
equations with cosmological constant, A 

Ef,u := + Ag^^ - |r^^ = 0, (A.38) 

are not written just in terms of the invariants ()A.34p . However, notice that for any symmetric tensor t^jy 

= + 2ny,)o + 2,5'^^,„(^^fe + 25"''ni(^t,),. (A.39) 
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Recalling (|A.35p we see that the tensor 

V := V + iB- E%, + 2{B - i?'),(M*.)o " 2<5'^^S - 25'^^F,,(^t,), (A.40) 

is gauge invariant. Now if we choose = E^i, then tf^u = 0, since the FRW metric solves the Einstein 
equations. Therefore E^i, = E^^, is gauge invariant, and in the second form it can be written just in terms 
of invariants when decomposing E^y = G^y + Ag^^ — kT = 0. The same holds for the perturbations. 
Explicitly, for instance, 



5Goo = SGoo + {B-E')G'oo + 2{B-E'yG 



00 



1 



SGoa = 6Goa + [B-E']^Goo-i:S"'{Fa + E,ayG, 



2 



cd 



SGab = 6Gab + {B-E')G,,--6'''[F^a + E,(a]^,^G,d, (A.41) 

where we used Goo = 3^2, Goa = 0, Gab = -(2W + n'^)6ab- Equation fOTI) equals (fXiOll with G^^ 
replaced by t^jy, whenever t^u is diagonal. Notice that the additional terms in (lA.4ip only contribute to 
scalar and vector perturbation contributions. Hence 6G^u = SG^u for tensor perturbations. We will see this 
explicitly. 

Applying ()A.40p to i^jjy = g^^ we find 

(5^00 = 2a2$, 5gQa = a^Va, Sgab = 2a^^5ab (A.42) 
while, applying it to t^^i, = T^y, we find after some algebra 

<5Too = ^(2C!Z' + a\'{0Z-2v{C)aH 



A 

Sfab = \ 



Xpa^hab + 2mab + ^Sab (2(C?$ + 2c' Z' - a\'{OZ 



(A.43) 



The computation for the Einstein tensor itself is more complicated and we divide it into modes. 



A. 5 Tensor perturbations 

For tensor perturbations we have 6 goo = 5goa = ^'^'^Sgac,d = ^'^'^^Qcd = and 5gab = a?hab- It immediately 
follows from (|A.22p that 5Roo = SRoa = and 

SRab = [^^9ab^ - 2W (y-^^Sab^ + ^^^'ab " ^^^9^b, (A.44) 

from which we immediately infer 6R = 0, so that 6{gabR) = R^gab- Using R = 6(W + TC^)/a'^, we find 
(5 Goo = SGoa = and for 5Gab = SRab — ^R^dab with 6 gab = o?hab after some algebra 

5Gab = \ [Kb + 2^/1',^ - t.hab) - {2n' + n^) hab- (A.45) 

As already mentioned, 5G^u = 5G^y for tensor perturbations. 

The term in the second square bracket of (|A.45p is not displayed in [9]. However, notice that there one 
perturbs the mixed components Gy = g^^^Gpy. Thus 

5Gii = -r^Sgaxt'Gpy + r'^Gpy . (A.46) 
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Hence <5G[! = 5^° = and 



5GI = -^^dgacGcb + ^5Gab = ^ [-habd'''^ + 5Gai}j = ^ {K, + 2nh',, - Ahab) , (A.47) 

where we used that d'^'^Gcd/'^ = — (2W + Ti."^). This is the equation we find in [S]. The reason why we 
display here the twice covariant tensor components is that it is this form that one finds more directly in the 
Hamiltonian formulation. 

A. 6 Vector perturbations 

In this case Sqqq = 0, dg^a = o?Sa, ^Qab = 2a^F(„ ;,'). Since both 5a, Fa are transversal, we have 6'^'^6goc,d = 
5'^'^6gcd = 0, so that we immediately find from (|A.22h that SRqq = 0. The remaining equations in (IA.22P 
simplify to 

SRoa = -^AVa+{n' + 2n^)Sa 

SRab = -(v(^^f,^ + 2nV^a,b))+2{n' + 2n')F(^a,b), (A.48) 

where we have used the linearly gauge invariant variable Va = Sa — F^. Next from ()A.23P we immediately 
see that SR = so that S{g^iyR) = R5g^y. Hence 5Gqq = and 

5Goa = -\AVa+(6^''^ya 

SGab = -(^(U)+2W(a,6))+2('5='^^)F(a,6) (A.49) 
where 5'"^Gcd/3 = -271' + U"^. Comparing with (TOTl) we see that 5Gqq = and 

5Gab = -Vl^^,^-2nVi^a,by (A.50) 

These are also the equations for the mixed components of the gauge invariant Einstein tensor that we find 
in [9]. 

A. 7 Scalar perturbations 

Now we have ^^oo = 2(j)o? 6goa = B^aO?, Sgab = 20^ {ip6ab + E^ab)- Hence d'^'^dgocd = a^AB, d'^'^Sgcd = 
2a^{^ip + AE). Thus (|A:22|) simplifies to 

5Roo = -A$ - 3 (V" + 7^ (V- + 0)') (A.51) 

5Roa = -2[^' + n^]^^ + m'{B- E')^a + {n' + 2n'^)E[^ 

5Rab = [-^ + <^> + (n' + 2n^)E] + (-A^ + V" + 2{n' + 2n^){'ip + ,/.) + 5W^' + n^') 6ab, 

where we again used the gauge invariant variables ^ = ip + Ti.{B — E') and $ = — T-i{B — E') — [B — E')' . 
As an intermediate result we have 

6'"^SRcd = A[-^ + <!> + {n' + 2n^)E]+3{-A^ + tP" + 2{n' + 2n^){i; + (l)) + 5ntp' + n^'). 

(A.52) 
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(A.53) 



5Gab = [-^ + ^-2{2n' + re)E]^^^ 

+ [A(^ - - 2{V + {2n' + n'^){'^ + + n{2^ + $)') + 2{n" + Hn'){B - e')] Sah 

where we used Gqo = SW^, (^^'^Gcd/S = -[2n' + n\ 
Hence, by (|X4T]) . we get 

^Goo = 5Goo + G'oo(i?-i^') + 2Goo(i3-ii;')' 
= 2{-/\^ + m^') 

SGoa = ^Goa + Gqq [B — E'] — -d^'^'GcdE'g^ 



involved, as one can explicitly check. 

B Linear Perturbations following an Alternative Route 

In the main text we computed the linear perturbations of the general equations motion for our manifestly 
gauge-invariant configuration observables, the metric Qjk and the scalar field H. That is, we used the second- 
order temporal derivative form of the equations of motion, in which the canonical momenta were eliminated 
by using the Hamiltonian equations of motion. This is rather tedious, because one cannot use the special 
symmetries of the background when eliminating the momenta, but has to assume a general background. 
As shown in the appendix of our companion paper [T], to linear order one can get those equations also 
by perturbing the general Hamiltonian equations of motion and then eliminating the perturbed momenta 
by using the perturbed equations of motion. This is shorter because one can use the properties of the 
background at an earlier stage in the computation. On the other hand, the procedure followed in the main 
text quickly becomes more economic for higher than linear order. 

In this section we will carry out this alternative derivation as a check of the result obtained in the main 
text. 

B.l Gauge invariant FRW equations 

We will need the gauge invariant FRW equations in Hamiltonian form: 

The metric takes the form Q^j. = A^5jk, where A is the observable scale factor. The corresponding extrinsic 
curvature is given by 



2 + 



A 





(B.l) 



The momentum P' conjugate to Qjj. is given by 



= ^det(Q) (g^"^Q'" - Q'^Q^'^) K, 



2AS^^ =: 16^''. 



(B.2) 
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The background scalar fields H, 11 are simply spatially homogeneous where 11 is given by 



The symplectic potential iq 



Thus, with 



G 



n = ^det(g)(: 



Ki A 



-A {61 A) + \eoU. 

K A 



J 



J ■= QIA = -12AA, I = — = -2A, 

QA 

we see that the non-vanishing background Poisson brackets are 

{J,A]- = K, {n,H}- = A. 



It follows that 
and 



-2DkP) + liD.E = 



kC. 



geo 



det(g) 

-- + 2^3 A 
2 

-— - + 2^3^ 
2AA 



xa 



matter 



1 




2A 


^det(g) 







+ Jdet{Q) (Q- 



'•k + v[ 



Ml 



(B.3) 

(B.4) 
(B.5) 

(B.6) 
(B.7) 



(Bi 



Here we have introduced the background matter energy density p. Notice that H = C = Cgco + C'matteD so 
that C is also the background Hamiltonian. 

The Hamiltonian background equations of motion are 



A 

ri 

j 



dC 



(B.9) 



J 



'dJ_ 
A— 

dC 



VIA 
1 



J2 



1(24^2 



3k. n 



J2 



(24^12 



where we have introduced the background matter pressure p. 



This involves spatial averaging over the general symplectic structure, or, alternatively, we only use the integrand. The 
justification for this is that the corresponding equations of motion derived via the Poisson brackets correctly reproduce the 
connection between A, H and the momentum conjugate to A, H displayed in (|B.2|I and (IB.3p . respectively. 
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The first two equations in (jB.OP correctly reproduce the equations 11 = A^'E and J = —12AA or / 
Taking the second time derivatives of A, H and using the last two equations in ()B.9p . we find 



-2A. 



A^ 

^1= 



A 

'a 



A 



P + 



e 

I3 



(B.IO) 



Here we have used the conservation law that the total energy density H = C =: e > is a constant of 
motion of the background Hamiltonian (it is not a constraint). The conservation equation C = I can be 
solved for A^ when using J = —12 A A and one finds 



(B.ll) 



which displays p^ust = —e/A^ as background dust energy density while Pjust = 0, since the dust is pressure- 
less. 

Equations (jB.lOp and (jB.lip are the familiar FRW equations. However, they now describe the physi- 
cal evolution with respect to the physical Hamiltonian H of observable quantities, rather than the gauge 
transformations of non-observables as is the case in the usual formalism. This happens due to the de- 
parametrisation through the dust whose only fingerprint is in the additional terms proportional to e. 



B.2 Linear perturbations of the Hamiltonian equations of motion 

We start from the general Hamiltonian equations of motion for our physical degrees of freedom which we 
display once again below: 



N 



-.U + C 



H 

Qjk 



Vd^t(Q) 
2N 



N 



Q, 



N 



.3 2 

Gjkmn P i.^]^Q)jk 
jm Tikn 



Vdet(Q) v'{E)+C^Il 



, , 2PJ™P'=" - pj^p^A + -Q^k C - TditQ Q^^ (2A + ( 

Vdet Q\ J 2 ^ 2\ 



+ Vd^[G-i 



-l^jkmnl (^D„,DnN) - N Rmn[Q]) + — iV^d^tQ H'^H'^ 



(B.12) 



with N = Y^l + Qi'^NjNk, Nj = -Cj/H, H = ^JC^ - Q^'^CjCk- Notice that in our convention C = Cg^o + 
Cmatter and Cgeo Contains the cosmological constant term. The bimetric Gjkmn = Qj{mQn)k~^/'^QjkQmn has 
the inverse (G"!)^'^"^" = [Qj{^Qn)k - QjkQmn], so that Gjfcpg(G-i)P'?™" = (5[^5^). In what follows, we again 
adopt the notation of the main text and denote background quantities by an overbar, while perturbations 
are denoted by the symbol 6. For instance, 6Qjk = Qjk — Qjk- 

In perturbing equations (IB.12|) around the FRW background we make first the following observations 
which will drastically simplify the subsequent analysis: 



1. In linear order the perturbation of the lapse vanishes 

1 



6N 



4= ( Q'^Cj6Ck + UQ^^Cfik 



0, 



(B.13) 



since C,- 



Hence Cj = due to spatial homogeneity (in dust space S). 
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2. The same holds for the second covariant derivatives of the lapse 



5D,DkN = S {djdkN - T]ldmN) = d,dk5N - I%dmSN + 6r]ldmN = 0, (B.14) 



since both = 0, DmN = due to spatial homogeneity. 

3. The terms quadratic in spatial derivatives of the scalar field have vanishing variation because at least 
one of the spatial derivatives is not varied and then vanishes. 

4. In the variation of the Ricci curvature terms only the term involving the variation of the Ricci curvature 
Rjk, which is linear in the Christoffel symbol, survives. 

5. As shown in the appendix of ^Lj, the linear perturbations of the conserved quantities H(a), Cj{a) and 
thus Nj = —Cj/H are conserved with respect to the linear equations of motion which are generated 
by the second order term of the perturbation of the physical Hamiltonian H. 

It follows that we can set everywhere N = N = 1. Notice however, that, while 6N = 0, Nj = 0, we have 
5Nj / 0. 

With these preparations out of the way, we can now perturb ()B.12p . In order to keep the formulae simple 
at intermediate steps we define SQjk = L5jk, 6P^^ = 15^^. Later we will substitute L = A?, I = —2A. 
Remember that P-'^, 11 are densities of weight one, so that (A^-' = Q^^Nj.) 

c^p^'' = d^{N"'p^'')-2N%p^^^, £^n = a„(iv™n). (b.is) 

Since A^, = 0, we have 



6Cj^E = 



21— 

'-j-Gjkmn i^^rnl^r^ 



L 

6/:^U = dn,{6N"'Il) = j[6N^]^^. (B.16) 



Here we have introduced the flat bimetric Gjkmn = ^j{m^n)k ~ ^f^^jk^mn with inverse Gj^^,^ = 5j(^^6n)k — 
SjkSmn, so that GjkpqGpqmn = 5j{m^n)k- Also we will use Einstein summation convention in what follows, ir- 
respective of index position on Kronecker (5's, that is, we define Q^^"^^ ■= Gjkmn and [G ^y^^"^ — [G ^]jkmn, 
as well as 6^^ = 6jk- However, notice that our fundamental perturbation quantities are Qjk,P-^^,Nj with 
that index position, and one has to take care of the additional metric contractions involved when the index 
is a priori not in that position. For instance. 



5Q 



dk 



Q Q 



We also often need 



(5Vdet(Q) = 

1 1 



i^det(Q)g"^" 



1 



7d^t(Q) 



(B.18) 
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Then we find after straightforward, but tedious calculations (simply using the product rule for linear varia- 
tions all the time) 
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6tl 



5U 



n 



L3/2 2L5/2 
1 



5Qr 



Li/2 _ 
-—v'{E)6Qr 



2 ^ 
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'~^jkmn 

"3/2 
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{™J,n) ^1/2 
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^jkmn'-'^ 



8 L 2 4^ 
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1 



5/2 ^ K ' 

-4T + ^-2^ 



(B.19) 



Here A = 6"^'^dmdn is the flat Laplacian and all derivatives appearing are flat derivatives. As in the previous 
subsection, we usedp=[n /L^ + v{E)]/{2X) andp= [n /L^ -v(E)]/{2\). The variation of the Ricci tensor 
is given explicitly by 

SRjk = ^ {'^SQm{j,k)m — ^^Qjk — SQmmJk) ■ (B.20) 

B.3 Second time derivative form of the perturbed equations of motion 

The explicit inversion of the first and third relation in (|B.19P for the perturbed momenta in terms of the 
velocities yields 



6U 



r, 



2LV2 



2\6N, 



(mj,n) 



^SQjk- 



(B.21) 



Taking the second derivative of the first and third relation in (jB.lOp . we substitute for SIl, 6P^^, using the 
second and fourth relation in ()B.19p . respectively. Afterwards we substitute for 511, 5P^^ , using (|B.2ip . In 
doing that, one has to remember that d/drSNj = 0, i.e. 5Nj is a constant of motion. 
We find again after some tedious algebr4§ 



3L 

-2L^' 



-v"5E + -A5E + 
2 L 
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2LV2 2L5/2 4L 



n 



2L5/2 



5Qr 



+ 



n 



[5N„ 



L 



1 ^ 



Li/2 



.L 



L 



(B.22) 



It is a good check to verify that the dimensionalities of the various terms match: in our convention, the 
spatial coordinates are chosen to be dimensionfree, while r has dimension of length. Thus L has dimension 
cm2 and A has dimension cm^. The scalar field E is dimensionfree, hence 11 oc A^E has dimension cm2. 
The potential term v{E) has the same dimension as iP/L^ which is cm~2. Likewise, Qjk has dimension cm2 
while Kjk oc Qjk/N has dimension cm-*^ so that P^^ oc det{Q)Q^^Q'^"'Kmn is actually dimensionless, just 
like / = -2A. 



'One should really write / {dT)'^5Qjk, rather than SQjt etc. However, these two quantities are numerically identical. 
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In the last line of ()B.22p we still must insert the last relation of (|B.19p . Since there are various bimetric 
contractions involved, we notice the identities 

GjkpqGpqmn — ^j{m^n)k '^^jk^mm Gjkmn^mn — "^^jk- (B.23) 

The calculation is very tedious but the result is rather simple. We directly substitute L = and / = —2A 
and find 



A ■ A 

SQjk = ^<5Qifc + 2— [(5iV(j]fe) + 2 



a) -I + (^-2^) 



^Qjk — '^^Rjk 



2 2 A 



We can simplify this expression further by making use of the background equations (|B.10p and (jB.lip which 
imply 

Therefore the SQmm term in (|B.24p vanishes and the 5Qjk term simplifies, leading to 

SQjk = jSQjk + 2^[5iV{i]fc) + 2j5Qjk - 26Rjk + 6jk Q^i?™™ - ^^'(^'^^ " (B.26) 

Equations ()B.22p and ()B.25P are precisely equations (j2.17p and (j2.2ip , derived by the more general formalism 
of section [2l 
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